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(57) Abstract 

The present invention provides a method for generating a 
representation of the extent of relatedness between at least two 
classes of cells. The invention also provides a method for generating 
a representation of the correlation between a first class of cells and a 
second class of cells. The correlation reflects a change in the nature 
and amount of nucleic acids present in the classes. In these methods, 
the cells in each class are chosen from among cells of a given cell 
type, cells from a given tissue, and cells from a given organ. The 
methods establish similarities or differences between the classes by 
defining a plurality of pairs of nucleotide subsequences, each pair 
consisting of a first subsequence and a second subsequence, and, 
in the nucleic acid of each class of cells, determining the presence 
of a fragment with the first subsequence at one end and the second 
subsequence at another end and having a length separated by the 
first and second subsequences, as well as a quantitation of the extent 
to which each fragment is present. The methods then determine the 
extent of relatedness reflecting the similarities or differences among 
the classes. The invention further provides display means displaying 
a representation of the extent of relatedness between the classes of 
cells, and displaying a representation of the correlation between 
the first class of cells and the second class of cells. Additionally, 
the invention provides a representation of the extent of relatedness 
between the classes of cells, and representation of the correlation 
between the first class of cells and the second class of cells. 
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GEOMETRICAL AND HIERARCHICAL CLASSIFICATION BASED ON 

GENE EXPRESSION 

FIELD OF THE INVENTION 

This invention relates to representations of the extent of relatedness between cells, cell lines, 
tissues, organs, or expressed sequences based on a genomic analysis of gene expression using software 
algorithm based analysis. 

RELATED APPLICATIONS 

This application claims priority to both United States Application Serial Number , 

filed September 16, 3999, entitled "GEOMETRICAL AND HIERARCHICAL CLASSIFICATION BASED ON 
GENE EXPRESSION", and United States Provisional Application Serial Number 60/101,009 filed 
September 17, 1998, entitled "PlIYLOGENOMICS AND PHARMACOGENOMICS", which are incorporated 
herein by reference in their entirety, 

BACKGROUND OF THE INVENTION 

The rapid development of genomics and proteomics in recent years has led to a burgeoning of 
applications making use of the new information provided. A significant area in which such information 
has been put to use is in the grouping and characterization of pathological states according to the 
differential expression of genes in such states. A corollary application is in grouping and characterizing 
the therapeutic effects of known or candidate pharmaceutical agents used in treating various pathologies. 
Algorithms employing a variety of statistical procedures have been employed to create heuristic displays 
of the information obtained from such analyses. These displays include large two dimensional, or even 
higher dimensional arrays in which the elements are coded, for example by false color coding, to 
represent a particular experimental result. Alternative displays include those in which the experimental 
data is used to generate cladistic or radiating tree structures as a representation of relatedness. 
Furthermore, it is also possible to use similar methods to group expressed sequences according to 
patterns of co-expression over several different biological states. 

For example, a system of cluster analysis for genome-wide expression in the yeast 
Saccharomyces cerevisiae and in primary human fibroblasts has been presented by Eisen et ah (Proc. 
Natl Acad, ScL USA 95:14863-14868 (1998)). In the yeast work, DNA microchip arrays carrying 
essentially every ORF from this organism were used. Differential expression was studied by varying the 
physiological state, including the diauxic shift, the mitotic cell division cycle, sporuiation, and 
temperature and reducing shocks. The human fibroblasts were stimulated with serum following serum 
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starvation, and examined using a microarray with 9,800 c DIM As representing approximately 8,600 
distinct human transcripts. Additionally, a further independent variable in these experiments is the time 
at which an assay point was taken. Data reflecting the differential gene expression in the various studies 
were analyzed using pairwise average-linkage cluster analysis (Sokal et aL, Univ. Kerns. ScL Bull. 
5 38:1409-1438 (1958)), which was used to compute a dendrogram that assembles all elements into a 
single tree 

Colon adenocarcinoma from 40 tumor samples were compared with 22 normal colon tissue 
samples using Affymetrix DIM A chips to which sequences from human cDNAs were bound (Alon et al, 
Proc. Natl Acad, ScL USA 96:6745-6750 (June 1999)). 3,200 full-length human cDNAs and 3,400 ESTs 

10 are represented in sets of 25-bp fragments, as well as such sequences containing a single base mismatch 
in the center of the sequence. The gene expression in both the tumor tissue samples and the normal 
colon samples, was assessed by hybridization. The statistical significance of the correlation between 
genes was assessed by calculating pairwise correlation coefficients. The clustering of the expressed 
genes was evaluated using an algorithm based on deterministic-annealing (Rose et ah, Phys. Rev. Lett. 

15 65:945-948 (1990); Rose, Proc. IEEE 96: 2210-2239 (1998)) to organize the data in a binary tree. Data 
are presented as a large two-dimensional color coded array, with genes displayed along one dimension 
and tissue samples along the other; artificial color values are assigned at each array point to indicate the 
extent of expression in a third dimension. Clustering analysis reveals patterns in the color distribution 
within the array which is disrupted when various randomization procedures are applied . The clustering 

20 of the genes in the data set reveals groups of genes whose expression is correlated across tissue types. 
The algorithm separated the tissues into distinct clusters. 

Pharmacological effects of compounds actually used or being screened for use in cancer 
chemotherapy were analyzed by cluster analysis at the National Cancer Institute (Weinstein et al, 
Science 275:343-349 (1997)). More than 60,000 compounds were screened against a panel of 60 human 

25 cancer cell lines. A 50% growth-inhibitory concentration of a compound in a given cell line, when 
analyzed across all cell lines, provided detailed information on mechanisms of drug action and drug 
resistance. Patterns of activity were first analyzed by the COMPARE algorithm (Paul! et al, J. Natl 
Cancer Inst 81 : 1 088 (1989); Jayaram, Biochem. Biophys. Res. Commun. 186:1600 (1992); Paull et al, 
In: CANCER Chemotherapeutic AGENTS, Foye (ed.), American Chemical Society, Washington DC, 

30 1993,pp. 1574-1581; Boyd et aL, Drug Dev. Res. 34:91 (1995)). The procedures developed rely on 

three databases, an S database characterizing structural information on the candidate compounds, an A 
database related to the 60 cell lines and a T database including information on molecular targets of 
action. In an example of the results of the analysis, a three dimensional array displaying compounds 
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versus targets, with a false color code providing a correlation coefficient in a third dimension for each 
position in the array, was developed. 

Certain problems arise upon consideration of the procedures currently in use for the correlation 
and clustering of genome- derived attributes. Use of DMA microchips inherently limits any analysis to 
5 the sampling of the DNA sequence fragments employed as the capture probes bound to the chips. 

Detection of any DNA fragment which does not hybridize with one of the capture probes is not possible, 
so that positive results are potentially lost. Additionally, a mutation or other allelic polymorphism may 
not bind to the capture probe under conditions of moderate or low stringency, so that again information 
relating to a positive result may be lost. 

10 For these reasons there is a need for methods of genomic statistical analysis based on more 

comprehensive accessibility to the genomes of the organisms being studied . Furthermore there remains a 
need for ways of presenting the information obtained in genomic analyses of relatedness of genes, and in 
genomic analysis of response to actual or candidate pharmaceutical agents, that includes information 
gleaned from a comprehensive access to the genomes in question. The present invention addresses these 

15 needs, for use is made in the invention of partial and full genomic sequences available from a large 
number of sequence databases in clustering analysis of the components appearing as independent 
variables in a particular study. 

SUMMARY OF THE INVENTION 

The invention provides novel methods of geometric and hierarchical classification between at 
20 least two classes of data sets. Data sets may represent cells, nucleic acid sequences, polypeptide 
sequences, or the like. The invention is able to utilize both standard DNA microchip arrays and 
non-DNA chip technology to provide input information on nucleic acid moieties of the specified classes 
of cells. The data are then treated in various ways to provide representations of relatedness that are 
readily interpretable by the human eye. The invention additionally provides novel methods for 
25 generating a representation of the correlation between at least two classes of cells, the correlation 
reflecting any changes in the composition and amount of nucleic acids present between the classes 

The cell classes may be from different sources for use in comparing differences between various 
cell populations. These differences include, but are not limited to, species differences, tissue differences, 
disease state differences, and drug treatment differences. Computer algorithms analyze input data 
30 reflecting differences between chosen cell classes and represent them in a meaningful way. 



WO 00/15851 -4- PCT/US99/21525 

Prior to the present invention, input information was obtained only using DNA-chip technology 
to analyze the nucleic acids of the cell classes to be compared. Drawbacks to these methods are that 
identifier sequences need to be already known and isolated, chip technology has size limitations related 
to the number of the nucleic acids immobilized on the chips, and, once the chips were manufactured, it is 
5 virtually impossible to expand nucleic acid parameters* The invention provides the use of 

GeneCalling™, a non-DNA chip technology, to assay differences between input cell classes . An 
unexpected result is that GeneCalling™ is able to provide sensitive comparisons between disparate 
groups above, thereby sidestepping the limitations inherent in the use of DNA chip technology when 
assaying input nucleic acid population. 

10 The invention provides a novel method for generating the extent of relatedness reflecting 

similarities or differences in the presence and quantitation of the fragments among the classes by 
calculating a distance that reflects the amplitude of a difference vector. In a significant embodiment of 
this method for generating the representation of relatedness, the extent of relatedness is provided by 
generating a tree structure reflecting the relatedness between any two classes „ The branches of the tree 

15 structure reflect the difference vectors and are ramified from nodes. 

The invention also provides a novel method for generating a representation of the correlation 
between classes of data sets. In a significant embodiment of the method for generating a representation 
of the correlation, the correlation is related to a set of orthonormal eigenvectors. In another significant 
embodiment of the method for generating a representation of the correlation, the representation is a 
20 cluster diagram or a dendrogram, and includes a tree structure reflecting the relatedness of the pathways 
involved in the biochemical or physiological response to a difference between cells of the two classes. 

The invention additionally relates to providing geometrical representations of differences 
between classes of data sets. The geometrical representations encompass, by way of nonlimiting 
example, principal component analysis and principal factor analysis, as well as reduced dimensional 
25 representations derived from them. The geometrical representations are based on differences determined 
between classes of cells using any method of analyzing for the presence of genes, nucleic acids, or 
fragments thereof, including nucleic acid microchip arrays and differential display of expressed genes or 
nucleic acid fragments. 

The invention also provides display means for displaying the representation of the extent of 
30 relatedness, the correlation, and the geometrical representations of differences between classes of data 
sets, as well as the representations themselves. 
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BRIEF DESCRIPTION OF THE DRAWING 

Figure 1 is a schematic flow diagram illustrating the principal steps involved in generating the 
various representations of the invention starting from a set of subsequence- selected fragments found for 
the samples. 

5 Figure 2 is a schematic flow diagram illustrating the primary steps involved in carrying out a 

principal component analysis. 

Figure 3 illustrates hierarchical clustering of four drugs with sterile water as an outgroinx 

Figure 4 is a graphical projection of drug treatments and controls onto principal factors. 

DETAILED DESCRIPTION 

1 0 The present invention relates to methods for preparing representations of the relatedness between 

cells of any two or more different classes of cells. The classes broadly encompass cells arising in animal 
and plant organisms, the cells further being normal cells or cells in a diseased state, including tumor 
cells. They further include cells that have been treated with a putative pharmaceutical agent The 
representations are obtained using experimental data that provide size and sequence information on 

15 nucleic acid fragments derived from each of the cellular sources. The fragments may be prepared from 
the nucleic acid content of the cells in each class in any of several ways. For example, in a particularly 
important embodiment, they may be subjected to digestion by particular pairs of restriction 
endonuc leases; alternatively, in another important embodiment, cell extracts may be subjected to 
amplification using specially designed primer oligonucleotides. The present invention also relates to 

20 methods for preparing representations of the relatedness in terms of co-expression between the nucleic 
acid fragments so produced. 

The invention further relates to the representations provided by these methods, and to display 
means on which such representations are displayed. The methods for preparing the fragments, such as 
the use of restriction endonucleases or the application of amplification primers, are chosen to provide 

25 subsequence information relating to the ends of the resulting fragments, while size determination 

provides the length of the fragment. In certain applications of these types of information, the size and 
subsequence results can optionally be scanned against databases providing known nucleic acid sequences 
in order to provide the identity of one or more candidate fragments of known complete nucleic acid 
sequences having the correct length and terminal subsequences (U. S. Patent No. 5,871,697; Shimkets er 

30 al 1999 Nature Biotechnology 1 7:798-803). This database look-up step is not a required feature of the 
current invention. For this reason, the present representations and methods are more comprehensive and 
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more informative of genomic variations among the samples than those currently known. As described in 
the Background of the Invention, currently known procedures are restricted in their comprehensiveness 
to those nucleic acid fragments that are applied to DNA microchips as probe sequences in a given 
procedure. Except for a narrowly limited set of model organisms with known genome sequence, the 
5 number of such probe sequences is considerably fewer than the number of known nucleic acid sequences 
available in sequence databases and employed in the present invention. Furthermore, even for fully 
sequenced genomes, genetic variation is not adequately probed with existing DNA microchips. This 
distinction characterizes an important advantage of the instant invention. 

The invention additionally relates to providing geometrical representations of differences 
10 between classes of cells- The geometrical representations encompass, by way of nonlimiting example, 
principal component analysis and principal factor analysis, as well as reduced dimensional 
representations derived from them. The geometrical representations are based on differences determined 
between classes of cells using any method of analyzing for the presence of genes, nucleic acids, or 
fragments thereof, including nucleic acid microchip arrays and differential display of expressed genes or 
15 nucleic acid fragments. 

As used herein, "sample" relates to a particular experimental state for which all the variables 
being studied in a project are held fixed. By way of nonlimiting example, if a variable is a class of cell, 
the "sample" refers to a particular cell type; if a variable is the subsequence pairs employed in the 
project, a "sample" refers to a particular subsequence pair; or if a variable is a set of putative 

20 pharmaceutical agents, a "sample" refers to a particular agent from the set.. As used herein, 

"representation" relates to any graphical, visual, or equivalent non-verbal display that provides an image 
of the results obtained according to the methods of the present invention. More specifically, a 
"representation" of the invention is obtained by transforming the quantitative results gathered by 
experiments underlying the invention. Examples of such data include, by way of non-limiting example, 

25 differential gene expression across classes of cell, and/or across a set of putative therapeutic agents, 
and/or equivalent types of experimental parameter. 

In important embodiments, a representation of the invention is generated by algorithms executed 
in a computer and is suitable for display on a display means, such as a display screen or monitor, 
employed in the operation of the computer. The representation is also suitable for storing in a storage 
30 module or data archive of such a computer. It is still further suitable for printing from the computer onto 
a medium such as paper or equivalent physical medium, and for recording it onto a portable storage 
medium, including, for example, magnetic media, CD ROMs and equivalent storage media. As used 
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herein, "display means" includes any of the objects and media identified above in this paragraph, as well 
as equivalent apparatuses and objects suitable for displaying the results of computational processes for 
visual inspection. 

As used herein, "extent of relatedness" is a characterization according to methods of the present 
5 invention of a degree of similarity or a degree of non-similarity between any two members of the same 
type of element; in particularly important embodiments, the type of element may be classes of cells. 

As used herein, a "putative pharmaceutical agent" relates to a chemical compound or a 
composition comprising at least one chemical compound which is a candidate for being a therapeutic 
agent Any such therapeutic agent may be used in treating a mammal suffering from a disease or a 

1 0 pathology, in treating the mammal with the therapeutic agent it is intended to attenuate the symptoms 
and/or the underlying causes of the disease or the pathology, to ameliorate the symptoms and/or the 
underlying causes, and/or to contribute to a cure of the disease or the pathology, Tsf on-limiting examples 
of a putative pharmaceutical agent include an agent drawn from a chemical compound library; an isolate 
from a natural source; a compound synthesized specifically as a putative agent; or a substance derived or 

15 obtained using the practices of genetic engineering and recombinant nucleic acid technology such as a 
recombinant protein, a fragment of a recombinant protein, a recombinant polypeptide, a fragment of a 
recombinant polypeptide, a recombinant peptide, or a nucleic acid including, for example an 
oligonucleotide intended as an antisense agent, and a recombinant gene intended for administration as a 
gene therapeutic agent. 

20 As used herein, a "fragment" of a nucleic acid relates to a contiguous portion originating from 

the genomic or cDNA-derived nucleic acid from a class of cells. The contiguous portion includes at or 
near each end a target subsequence defined according to the operational procedures disclosed herein, and 
includes all nucleotides in the sequence of the fragment bounded by the two target subsequences. The 
nucleotides between the two target subsequences, together with the subsequences themselves, define a 

25 "length" of the fragment, as used herein. The target subsequences are identified, for example, by 

contacting the nucleic acid from the cells with a specific pair of restriction endonucleases, or with a 
specific pair of oligonucleotide primers, and in equivalent ways. 

The information used in the present invention is obtained from experiments providing the results 
of differential gene expression wherein the difference relates to an experimental state and a reference 
30 state. Commonly a reference state refers to a normal, or an unperturbed, or a non -pathological class of 
cells. An experimental state may relate to a certain set of conditions applied to one class of cells, and the 
corresponding reference state then relates to the same set of conditions applied to a second class of cells.. 
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An experimental state may also relate to a class of cells in the presence of one or more putative 
therapeutic agents, in which case the reference state relates to the same class of cells in the absence of 
any putative therapeutic agent. An experimental state may furthermore be obtained from a class of cells 
that is of interest in a particular set of circumstances. This includes celts of a given cell type, cells from 
5 a given tissue, and cells from a given organ, and further includes cells that may be noncancerous or 
cancerous. Types of cell encompassed within the present invention include, by way of non-limiting 
example, endothelial cells, mesotheiial cells, and epithelial cells. Tissues and organs included within the 
present invention may be, by way of non-limiting example, lung, heart, skeletal muscle, smooth muscle, 
brain, central nervous system, peripheral nervous system, stomach, liver, kidney, reproductive tissues 
10 and organs, skin, and bone. Cancerous cells include, by way of non-limiting example, cells from prostate 
cancer, breast cancer, colon cancer, lung cancer, lymphatic or hematopoietic cancers, and also include 
cells obtained from tissue biopsies or from cell lines in the National Cancer Institute human tumor cell 
line panel. The cells subjected to analysis in the present invention may also originate from plants, yeast, 
fungi, and other taxonomic groupings. 

1 5 The methods of evaluating the extent of relatedness between classes of cells, for example, 

between a first class of cells and a second class of cells, are founded on evaluating the extent of 
relatedness of the expression of particular genes between the cells of the two classes. In a preferred 
embodiment of the invention, similarities and differences in the susceptibility of the nucleic acid present 
in the cells to digestion by specific pairs of restriction endonucleases are determined, according to the 

20 methods of the present invention, by procedures that are disclosed in detail in co-owned LL S. Patent No. 
5,871,697 to Rothberg et aL, and in Shimkets et al 1999 (Nature Biotechnology 1 7:798-803), both of 
which are incorporated herein by reference in their entirety. 

Briefly, for any experimental state of a class of cells, the nucleic acid content of the cells, 
preferably in the form of a preparation of cDNA from the cells, is subjected to restriction endonuclcase 

25 ("RE") digestion by specific pairs of endonucleases. Each member of the RE pair is chosen to optimize 
the likelihood that a restriction fragment resulting from the nuclease digestion will be a unique fragment, 
in an important implementation of this method, the restriction nuclease digestion is carried out on cDNA 
prepared from the cells of the class in the given experimental state. This implementation leads to 
emphasis on genes that are expressed in the experimental state, many of which may be characteristic of 

30 the given experimental state and be more poorly expressed, or not expressed at all significantly, in a 
different experimental state. A large number of specific pairs of nucleases may be employed. 
Alternatively, expression of a gene may be repressed in a characteristic way in a given experimental state 
and be expressed at a higher level, such as at a constitutive level, in a different experimental state* By 
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way of non-limiting example, several pairs of restriction nucleases that may be employed in 
implementing the present invention are disclosed in U. S. Patent No. 5,871,697. 

In an alternative embodiment, the extent of relatedness may be obtained by amplification 
fragment length polymorphism analysis ("AFLP"). Briefly, amplification of the nucleic acid content of 
5 the class of cells being examined is subjected to a primer-dependent amplification procedure in which 
any of a set of primer pairs is used to initiate amplification. Amplification procedures are described in 
considerable detail in, for example, Innis et al. 9 PCR PROTOCOLS, A GUIDE to METHODS AND 
Applications, Academic Press, New York (1989), and Innis et al, PCR STRATEGIES, Academic Press, 
New York (1995). The primers of each primer pair are different from each other, and reflect different 

10 subsequences that are the object of the amplification process. Amplification may proceed by any 

procedure, including polymerase chain reaction, known in the field of molecular biology.. In AFLP, the 
length of an am pi icon found in a given experimental state differs from the length found in a different 
experimental state. This may arise, for example, if the given experimental state arises from a mutation 
that occurs in a subsequence recognized by a primer used in the amplification reaction. It may also arise 

1 5 from a deletion from, or an insertion into, the nucleic acid of the cells in that state. 

The experimental and computational procedures that may be employed to generate the 
representations of the present invention are described generally below. 

Measurements 

At the outset, the gene expression levels are determined experimentally. This can be done, in a 
20 preferred embodiment, by following the general protocols of differential expression using restriction 
endonucleases (U.S.. Patent No. 5,871,697). For each pair of restriction enzymes and each biological 
sample, a pool of fluorescently- labeled DNA fragments is generated. Electrophoresis is then performed 
to separate these fragments based on size, and an intensity, designated as I srt (x), where s labels the 
sample, i.e., the cell class; r labels the restriction enzyme pair, i.e., the gene fragment; t labels the trial, 
25 and x is the length of the fragment as determined by electrophoresis, is detected. The length x may be 
either a continuous index or a convenient discretization. As an example, the resolution of the 
electropherogram may be set to a discretization of 0.1 nucleotide ("nt"). Commonly three independent 
trials are performed. A mean signal I sr (x) is then obtained by averaging over the n t trials, 

I sr (x) = (l/n I )Z t I srt (x)S i (I) 
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Next, lengths x for each restriction enzyme pair r where some of the samples have a significant 
difference in measured intensity are identified. Such a difference is determined with respect to cell 
types, or with respect to the presence vs. the absence of a putative pharmaceutical agent. Labeling the d th 
such difference d, the values l sd = I sr (x ) are then collected* Any of several methods for identifying 
significant differences may be employed, some of which are outlined herein. For example, an important 
method involves the following computational steps: 

1 . The mean l r (x) - L s I sr (x) is evaluated. 

2. All positions, i.e. lengths, where, for at least one sample, l sr (x) - I r (x) is larger than some 
threshold value, are marked. 

3. The largest value of I sr (x) - I r (x), determined as a difference between a sample state and 
the mean for restriction enzyme pair r, is found and the length x, indexing the difference, 
is marked. 

4. Step 3 is repeated for succeedingly smaller values of the intensity difference. If the 
length x that marks the current largest difference is within a distance w from the length 
of a previously identified difference, the current difference is skipped and the next 
smaller difference is considered. 

5. Step 4 is repeated until there are no more differences to consider. 

Another method involves finding differences that meet a statistical criterion. A particular 
example of such a method involves the computational steps of: 

1 . defining a set of sample classes and assigning each sample to a particular class c; 

2.. for each restriction enzyme pair r and length x, evaluating the F-statistic for the set of 
measurements I sr (x) and the classes c to which samples are assigned, thereby providing 
the probability p r (x) that any differences between sample classes may be explained by 
random variation (See, for example, P, Hinton, Statistics Explained, Routledge 1995); 

3. ordering the probabilities p r (x) from smallest (most significant) to largest (least 
significant); 

4. optionally truncating the list at some threshold value of p r (x) above which differences are 
no longer considered significant (accepted values are p r (x) - 0„01 to 0.05); 

5. finding the smallest value of p r (x) and marking the length x as a difference for restriction 
enzyme pair r; 
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6„ repeating step 4 and determining whether the length x that marks the current difference 
is in a region that is within a distance w of a previous difference, in which case the 
current difference is skipped and the next smaller distance is considered; and 

7. continuing until there are no more differences to consider. 

These exemplary computational procedures provide a set of measures of intensity l sd for the class 
of cells in sample s at difference d. 

Distances 

For hierarchical clustering, a distance D ss > may be defined as the distance in vector space 
between pairs of samples s and s\ A variety of methods for calculating D Sfi are available. Some 
examples, which are intended as being noniimiting, are provided below. 

D £S ' as a scaled correlation function: 

1. One calculates ju d = (l/n s ) X s I sd and a d = [( l/n s )S s (I sd - \i d ff 5 . If data is missing, for 
example no measurement of I sd exists for some sample s, that sample is excluded from the 
sum and n s is reduced by 1 . 

2. One calculates J S[ , = (l sd - |i d ) / ar d . If data is missing for I sd , then J sd is defined as J sd = 0. 

3. One calculates ii s - (I/n d ) S d J, d and <y s = [(l/n d )S d (J sd - txjT 5 . 

4. One calculates K sd = (J sd - ii s ) / a s . 

5 . One calculates the covariance matrix S ss = ( 1 /n d ) Z d K sd K Vd 

6. One calculates the correlation matrix C ss = S ss > / [ S ss S sV ] 0 ' 5 . 

7. One calculates D ss - [2 - 2 C ss . ]°" . 



D ss > as a Euclidean distance: D ss - [ Z d (I sd - l 5 . d ) 2 ] 05 . 



D ss as a Pearson distance: D ss - [ E d (I sd - I sM ) : / rj d 2 f 5 where a d is defined in step 1 of scaled 
correlation function above. 
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D ss - as a pairwise Pearson distance: 

1 . One calculates the covariance matrix S ss , = ( l/n d )[ S d I sd I 5 d - (E a I sd ) (L d I s . d ) / n d ]. 

2. One calculates the correlation matrix C ss . = S ss / [ S ss S ss ]° 5 . 

3 . One calculates D ss . = [ 2 - 2 C ss > ]° 5 . 



as a Mahalanobis distance: 

1 . One calculates the covariance matrix S dd - - (Z s I sd I sd ) - (S 5 I sd ) (L s I S(r )/n s 

2. One calculates the correlation matrix C dd . = S dd / [ S dd S d - d . ]° 5 and its matrix inverse C' l M 
3 , One calculates D ss * = [ Z d(r (I sd - I, d ) C l d(V (I S(1 - - W) ] 0 5 

It is contemplated that other distance methods known in the art may be used in the invention, 
such as Spearman correlation, and the like. Other methods known in the art can be found, for example 
and not be means of limitation, in K. V. Mardia, J. T. Kent, and J. M. Bibby, MULTIVARIATE ANALYSIS, 
Academic Press, New York, 1979. 

Hierarchical Clustering 

The distances can be used to perform hierarchical clustering of the samples. A genera! algorithm 
for clustering is described below. 

1 . Each sample s is assigned its own initial cluster c. 

2. One calculates all the distances between pairs of clusters and finds the smallest distance. 
These two clusters are joined into a single cluster and the number of clusters is decreased 
by 1. 

3. Step 2 is repeated until only a single cluster remains„ 

In order to implement this algorithm, a method to calculate the distance between pairs of clusters 
is also required. Some nonlimiting examples of such calculations, using well-known methods, are 
indicated below. 

Nearest neighbor, single linkage: The distance between clusters c and c' is the smallest distance 
D ss , where s ranges over all samples in cluster c and s 5 ranges over all samples in cluster c\ 

Unweighted pair group method using arithmetic averages (UPGMA), also known as average 
linkage: The distance between clusters c and c' is (2^ D ss ) / (n c n c ) where s ranges over all samples in 
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cluster c, s' ranges over all samples in cluster c\ n c is the number of samples in cluster c, and n c is the 
number of samples in cluster c\ 

Furthest neighbor, complete linkage: The distance between clusters c and c' is the largest 
distance D ss - where s ranges over all samples in cluster c and s' ranges over all samples in cluster c\ 

5 Other distance- based hierarchical clustering methods are well-known. See, for example, Wen- 

Hsiung Li, MOLECULAR EVOLUTION, Sinauer Assoc, 1997. 

Software packages are available to perform the clustering and display the results. See, for 
example, Phylip, Joe Felsenstein, http://evolution.genetics.washington.edu for clustering, and Treeview, 
Rod Page, http://taxonomv.zoologv.gla.ac.uk/rod/treeview.html for display. The source code for the unit 
1 0 within Phylip employed for the clustering, and the downloaded executable file of Treeview for Windows 
95 and Windows NT, as well as a manual for Treeview, are available from the owner of the present 
application. 

Two-Dimcnsional Clustering 

It is also possible to cluster the distances, rather than clustering the samples. One simply 
1 5 exchanges the roles of the samples and differences in the equations above. Furthermore, it is possible to 
perform clustering of both samples and differences, and then to display the measurements I sd in which 
both samples and differences are presented in cluster order. 

Principal Component Analysis and Principal Factor Analysis 

Principal component analysis is described in standard texts. See, for example, Mardia, Kent, and 
20 Bibby. To perform principal component analysis, one begins with a correlation matrix C ss * as defined 

above, in the section "Distances". (Alternatively, one could use the covariance matrix S ss .)- Eigenvalues 
and eigenvectors, defined such that C ss g s i = a^ , where the i th eigenvalue is ^ and its eigenvector is g si , 
are calculated. The eigenvalues are ordered from largest to smallest: ai > a 2 > ... > a s . To obtain a 
reduced dimensional depiction of the samples, a number of desired dimensions k is chosen. Then, in k- 
25 dimensional space, sample s is represented as the point (g st , g s2 , . . g sk ). Samples that are close in the k- 
dimensional space have similar expression profiles and may be considered to be related. 

As an alternative to using the correlation matrix C ss - as the starting point for principal component 
analysis, it is possible to calculate principal components using the inner product matrix from 
multidimensional scaling defined as 

30 B=HCH (2) 
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where C is the correlation matrix, H is the centering matrix with diagonal elements given by 1 - (1/n) and 
off-diagonal elements - (1/n), where n is the number of items being correlated. (See, for example, 
Mardia, Kent, and Bibby, Multivariate Analysis, and Arkin, Shen, and Ross, Science 277: 1275 (1997)). 
The k th principal component is then the k ,h eigenvector of B normalized to unit length and ordered by 
5 decreasing eigenvalue A, k , and the k ,h principal factor is obtained by scaling the eigenvector by \ k m . The 
projection of sample s onto the k th principal factor is the element of the factor for row s. The 
components or factors are ordered from 1 (corresponding to the most informative) to n (corresponding to 
the least informative)* By using some, but not all, of the components or factors, the samples can be 
represented in a small-dimensional geometric space. Furthermore, the amount of information retained in 
1 0 the representation can be related to the eigenvalues of the components that are used (See Mardia, Kent, 
and Bibby). 

A centered inner product matrix B appropriate for principal component or prinicpal factor 
analysis can also be obtained from any distance matrix D ss> as 

B = HAH (3) 

1 5 where 

A s , =-l/2(D ss 0 2 - (4) 

To perform principal factor analysis, factor i is defined as h si = a ; 0 5 g si where, as before, a; is the 
eigenvalue of the i th eigenvector g si . An orthonormal rotation matrix G (Ej GjjG kj is 1 if i = k and 0 
otherwise, det(G) - +1) is introduced and the factors are rotated to obtain rotated coordinates for the 
20 samples. Thus, to obtain a k-dimensional representation of the locations of the samples, the following 
operations are performed: 

1 . One calculates the correlation matrix C ss * or the co variance matrix S ss - where s and s' 
label individual samples. 

2. One calculates the eigenvalues a^ and eigenvectors g si for the matrix, with & } > a 2 >. .„> a s . 

25 3. Unrotatcd factor loadings h si = a^ 5 g sj are defined. 

4. The first k factor loadings and an orthonormal rotation matrix G are selected. The j th 
coordinate of sample s in the rotated space is Ej. h sJ - G rj 

The rotation matrix G may be optimized according to standard criteria. See, for example, 
Mardia, Kent, and Bibby, Ch„ 9.6 on Varimax rotation, supra. The rotated axes represent factors that 
30 influence the observed measurements for the samples. 
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In implementing the methods of the present invention, these operations may be sequentially 
combined in any of several ways according to the intended display, i.e., the nature of the relatedness that 
is intended to be shown. 

Also, the information from the principal factors can be used to help filter the experimental noise 
5 from the correlation functions. For example, it is possible to select a cut-off principal factor j < n, then 
compute distances and correlations between samples based on their representation in the j~dimensionaI 
principal factor space. 

As a nonlimiting example of the computational procedures that may be employed in the present 
invention, a schematic overview of procedures that may be adopted is presented in Figure 1 . The 

10 experimental results represent the sample-dependent and selection-dependent intensities obtained in an 
experiment, arrayed in a measurement matrix. In the implementation shown in Figure 1, the difference 
bands having various, defined, nucleotide lengths are arrayed as the columns of the matrix; they are 
obtained in various experiments that are selected using different members of the sets of subsequence 
pairs. The samples represent the classes of cells, or cells treated with a set of putative pharmaceutical 

1 5 agents, or analogous sample sets, and are arrayed as the rows. 

The values arrayed in the measurement matrix may then be subjected to correlation analysis to 
provide either direct sample correlations or correlations of differences. The measurement matrix can 
also be subjected to a calculation providing a vectoral distance between samples; such a sample distance 
may also be obtained from the sample correlation result. The distance vector can further be subjected to 
20 a linkage analysis to provide hierarchical clustering of the samples. Additionally, the correlated samples 
may be subjected to principal component analysis providing the principal factors contributing to a state 
or to a difference . 

A nonlimiting example of the way in which a principal component analysis may be carried out, 
using methods described herein, is presented in Figure 2. The correlation matrix or the centered inner 
25 product matrix described above is subjected to appropriate operations to provide the principal 

components and the principal factors, based on their eigenvalues and eigenvectors. Advantageously a 
reduction in the number of dimensions employed in the number of eigenstates may provide a filtering 
effect, reducing the noise in the vector distances calculated. 

The representations provided in the present invention find use in various applications of 
30 genomics in the biological and medical fields. Extents of relatedness and correlations provide rapid 
overviews of enzymatic reactions, metabolic pathways, and physiological effects that become 
distinguished when comparing states. When a pathological state is compared with a normal state, for 
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example in a mammal, and especially in a human, the display of distinguished pathways is instructive in 
the development of therapeutic approaches and/or therapeutic agents for the treatment of the pathological 
state. When a putative pharmaceutical agent is compared to a state that omits the agent, or when one 
such agent is compared with another, important information is provided relating to the metabolic 
5 reactions induced by or undergone by the agent or agents, leading to optimal choice of such agents. This 
information may also provide leads to the development of novel pharmaceutical agents. If the genome 
being studied is a plant genome, such as the genome of an important crop plant, analogous principles 
apply. 

Nucleic acid assays 

10 The present invention provides a method for generating a representation of the extent of 

relatedness between at least two classes of cells. In this method, the cells in each class are chosen from 
among cells of a given cell type, cells from a given tissue, and cells from a given organ. Generation of 
nucleic acids from the cell samples of choice may be as described in the GcneCalling™ methodology. 
See U.S. Patent No. 5,871,697. The method includes the steps of: (a) defining a plurality of pairs of 

1 5 nucleotide subsequences, each pair consisting of a first subsequence and a second subsequence; 
(b) isolating the nucleic acid of each class of cells and assaying for the presence of a nucleic acid 
fragment with the first subsequence at one end and the second subsequence at another end and having a 
length separated by the first and second subsequences, and quant itating the extent to which each 
fragment is present; and (c) determining the extent of relatedness reflecting similarities or differences in 

20 the presence and quantitation of the fragments among the classes using software algorithm programs 
known in the art. 

One important embodiment of this method, i.e., determining the presence of the fragments and 
quantitating the amounts present, as described in step (b) above, is carried out by a process that includes 
the steps as follow. First, samples of the nucleic acid from the cells of each class are digested with a 
25 plurality of specific pairs of restriction endonucleases ("REs"). Each sample is treated by one RE pair, 
where one RE of the pair targets the first subsequence described in step (a) above, and the second RE of 
the pair targets the second subsequence, with each digestion providing specific restriction fragments. 

Second, double stranded adapter DNA molecules are hybridized to the fragments. Each adapter 
DNA molecule comprises: (?) a shorter strand, preferably having no 5' terminal phosphate, consisting of 
30 a first and second portion, the first portion being a region at the 5' end that is complementary to the 

overhang produced by one of the REs of the given pair and a second portion hybridizable to the opposite 
longer strand of the adaptor, and (//) a longer strand, preferably having no 5' terminal phosphate, 
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consisting of a first portion at its 3' end complementary to the above-mentioned second portion of the 
shorter strand, and an optional second portion at its 5' end comprising a unique region not hybridizable 
to any sequence present in the original sample population. See U.S. Patent No. 5,871,697. The longer 
strand is optionally labeled with fiuorochrome 208, although any DNA labeling system that preferably 
5 allows multiple labels to be simultaneously distinguished is usable in this invention. See, e.g., Ausubel, 
et al Current Protocols in Molecular Biology, John Wiley & Sons, New York, NY, 1993. 

Third, output signals from each ligated fragment are detected for each sample population so 
treated. Each ligated fragment generates output signals that characterize (a) the presence of the given 
subsequences corresponding to the RE pair used in a particular run, (b) the length between the two 
10 subsequences corresponding to the two REs employed in a given run, and (c) the quantitation of the 
relative amounts present of each fragment so generated in a given run . 

Optionally, a nucleotide sequence database may be searched for sequences that are predicted to 
produce, or alternatively, not produce, the one or more output signals generated by the nucleic acid from 
the cells of each class, given the parameters described above. The analysis methods comprise, first, 

15 selecting a database of DNA sequences representative of the DNA sample to be analyzed, second, using 
this database and a description of the experiment to derive the pattern of simulated signals that would be 
generated, contained in a database of simulated signals, that will be produced by DNA fragments 
generated in the experiment, and third, for any particular detected signal, using the pattern or database of 
simulated signals to predict the sequences in the original sample likely to cause this signal. Further 

20 analysis methods present an easy to use user interface and permit determination of the sequences actually 
causing a signal in cases where the signal may arise from multiple sequences, and perform statistical 
correlations to quickly determine signals of interest in multiple samples. A sequence from a searched 
database is predicted to produce the one or more output signals when that sequence has both (a) the same 
length between occurrences of target nucleotide subsequences as is represented by the one or more 

25 output signals, and (b) the same target nucleotide sub-sequences that are represented by said one or more 
output signals, or target nucleotide subsequences that are members of the same sets of target nucleotide 
sub-sequences represented by the one or more output signals. 

A first analysis method is selecting a database of DNA sequences representative of the sample to 
be analyzed. In the preferred use of this invention, the DNA sequences to be analyzed will be derived 
30 from a tissue sample, typically a human sample examined for diagnostic or research purposes. In this 
use, database selection begins with one or more publicly available databases which comprehensively 
record all observed DNA sequences. Such databases are GenBank from the National Center for 
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Biotechnology Information (Bethesda, Md„), the EMBL Data Library at the European Bioinformatics 
Institute (Hinxton Hall, UK) and databases from the National Center for Genome Research (Santa Fe, 
N.Mex.). However, as any sample of a plurality of DNA sequences of any provenance can be analyzed 
by the methods of this invention, any database containing entries for the sequences likely to be present in 
5 such a sample to be analyzed is usable in the further steps of the computer methods. 

A second analysis method uses the previously selected database of sequences likely to be present 
in a sample and a description of an intended experiment to derive a pattern of the signals which will be 
produced by DNA fragments generated in the experiment. This pattern can be stored in a computer 
implementation in any convenient manner. In the following, without limitation, it is described as being 
1 0 stored as a table of information. This table may be stored as individual records or by using a database 
system, such as any conventionally available relational database. Alternatively, the pattern may simply 
be stored as the image of the in-memory structures which represent the pattern, 

A second important embodiment of this method, i.e., determining the presence of the fragments 
and their quantitation, as described in step (b) above, is carried out by a process that includes the steps as 

1 5 follow. First, for each pair of nucleotide subsequences selected, a pair of oligonucleotide primers are 
provides, the pair consisting of a first primer and a second primer, wherein the first primer is 
complementary to the first subsequence and the second primer is complementary to the second 
subsequence. Second, the nucleotide sequence between the first subsequence and the second 
subsequence are amplified using the oligonucleotide primers to prime the amplification, thereby 

20 providing an amplicon characterized by the subsequence pair, a length between the two subsequences 
corresponding to the two pr imers employed in each pair and a quantitation of the extent to which each 
amplicon is present. Third, output signals are generated as above for each amplicon, each output signal 
characterizing (a) the subsequences of the pairs of primers, (b) the length, and (c) the quantitation 
Optionally, a nucleotide sequence database may be searched for sequences that are predicted to produce, 

25 or alternatively, not produce, the one or more output signals generated by the nucleic acid from the cells 
of each class, given the parameters described above. Analysis methods are as described above. 

This invention can be applied, for example and not by way of limitation, to in vitro cell 
populations or cell lines, to in vivo animal models of disease or other processes, to human samples, to 
purified cell populations perhaps drawn from actual wild-type occurrences, and to tissue samples 
30 containing mixed cell populations. The cell or tissue sources can advantageously be a plant, a single 
celled animal, a multicellular animal, a bacterium, a virus, a fungus, or a yeast, etc. The animal can 
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advantageously be laboratory animals used in research, such as mice engineered or bred to have certain 
genomes or disease conditions or tendencies. 

Cells used in the invention may be obtained from a mammal, preferably a human, having or 
suspected of having a diseased condition. In one embodiment, the diseased condition is a malignancy, 
5 The in vitro cell populations or cell lines can be exposed to various exogenous factors to determine the 
effect of such factors on gene expression . In a preferred embodiment, the exogenous factor is a putative 
pharmaceutical agent. Cells so contacted with a putative pharmaceutical agent are treated with an 
amount of the agent sufficient to effect a change in the state of those cells or with an amount of the agent 
less than or equal to a predetermined upper limit of dosing concentration, prior to their being assayed, 
1 0 Measures of relatedness and extent of correlation may be made between cells so contacted with putative 
pharmaceutical agent and, for example, cells not so contacted. 

Extent of relatedness methodology 

The present invention provides a representation of the extent of relatedness between a first class 
of cells and a second class of cells. The cells in each class are chosen from among cells of a given cell 

15 typ^ cells from a given tissue, and cells from a given organ, as described above. The extent of 

relatedness reflects similarities or differences in the presence of pairs of nucleotide subsequences, each 
pair consisting of a first subsequence and a second subsequence, in a nucleotide length separating the 
first and second subsequences of the pair and in a quantitation of the extent to which each pair having the 
determined length is in the classes of cells. Input information of the fragments to be analyzed are 

20 obtained by methods of nucleic acid analysis and quantitation as described in the NUCLEIC ACID ASSAYS 
section above. 

The measure of relatedness is provided by calculating a distance that reflects the amplitude of a 
difference vector. A difference vector is defined as a difference between a first vector and a second 
vector. Herein, the first vector reflects information derived from the quantitation for each subsequence 
25 pair obtained for the first class of cells, and correspondingly, the second vector reflects the analogous 

information derived from the second class. The different elements of each vector relate to data obtained 
using different subsequence pairs. 

In an embodiment of the representation, the extent of relatedness is related to a distance. This 
distance reflects the amplitude of a difference vector that is a difference between a first vector which 
30 reflects information derived from the quantitation for each subsequence pair obtained for the first class 
and a second vector which reflects the corresponding information obtained for the second class. The 
different elements of each vector relate to data obtained using different subsequence pairs. 
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In an additional significant embodiment, the representation includes a tree structure reflecting 
the extent of relatedness is provided by generating a tree structure reflecting the relatedness between any 
two classes. The branches of the tree structure reflect the difference vectors and are ramified from 
nodes. 

5 In important embodiments of the representation of the extent of relatedness, the representation is 

obtained employing the methods of the invention, including the methods that have been summarized in 
the paragraphs immediately above. 

In additional significant embodiments of the representation of the extent of relatedness, the cells 
in at least one class are obtained as described in the NUCLEIC ACID ANALYSIS section above. 

10 Correlation analysis methodology 

The invention also provides a method for generating a representation of the correlation between 
a first class of cells and a second class of cells. The correlation reflects a change in the nature and 
amount of nucleic acids present in the classes. In this method, the cells in each class are chosen from 
among cells of a given cell type, cells from a given tissue, and cells from a given organ. The method of 
1 5 nucleic acid analysis and quantitation are as describe in the NUCLEIC ACID ASSAYS section above. 

Upon generation of a signal output, the correlation between the cells of the first class and cells of 
the second class are correlated, and a representation of the correlation is prepared. The quantitation of 
the fragments in the invention corresponding to the RE pair used in a given run and the length of each 
fragment so generated; thereby providing a quantitative measure of the extent to which the nucleic acid 
20 present in the cells in each class contains fragments having the specific subsequence pairs and the 
nucleotide length between the pairs. 

In a significant embodiment of the method for generating a representation of the correlation, the 
correlation is related to a set of orthonormal eigenvectors, as described in the DISTANCES section above. 
The elements of the basis set upon which the eigenvectors are constructed reflect particular biochemical 

25 or physiological pathways correlated between the cells of the two classes Each of these eigenvectors is 
associated with an eigenvalue that is an integer greater than zero. After defining an upper limit of the 
eigenvalues to be used, the coefficients of the basis set elements in each eigenvector whose eigenvalue is 
less than or equal to this upper limit reflects the contribution of the corresponding pathway to the 
biochemical or physiological differences correlated between the cells of the first class and the cells of the 

30 second class. 
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In another significant embodiment of the method for generating a representation of the 
correlation, the representation is a cluster diagram or a dendrogram, and includes a tree structure 
reflecting the relatedness of the pathways involved in the biochemical or physiological response to a 
difference between cells of the two classes. In obtaining this representation, a correlation matrix is 
5 calculated that provides a distance determination in which the distance reflects the amplitude of a 

difference vector. This vector is a difference between two vectors each of which reflects information 
obtained for the response of one of the two classes to the difference between the classes, and wherein the 
branches of the tree structure reflect the difference vectors and the branches are ramified from nodes. 

In additional significant embodiments of the representation of the extent of correlation, the cells 
10 in at least one class obtained as described in the NUCLEIC ACID ANALYSIS section above. 

Display means 

The present invention also provides a display means displaying a representation of the extent of 
relatedness between a first class of cells and a second class of cells. The cells in each class are chosen 
from among cells of a given ceil type, cells from a given tissue, and cells from a given organ, as 
15 described above. The extent of relatedness reflects similarities or differences in the presence of pairs of 
nucleotide subsequences, each pair consisting of a first subsequence and a second subsequence, in a 
nucleotide length separating the first and second subsequences of the pair and in a quantitation of the 
extent to which each pair having the determined length is in the classes of cells. 

In a significant embodiment of the display means, the extent of relatedness is related to a 
20 distance. This distance reflects the amplitude of a difference vector that is a difference between a first 
vector which reflects information derived from the quantitation for each subsequence pair obtained for 
the first class and a second vector which reflects the corresponding information obtained for the second 
class. The different elements of each vector relate to data obtained using different subsequence pairs. 

In an additional significant embodiment of the display means, the representation includes a tree 
25 structure reflecting the relatedness between any two classes, in which the branches of the tree structure 
reflect the difference vectors and the branches are ramified from nodes. 

In important embodiments of the display means displaying a representation of the extent of 
relatedness, the representation is obtained employing the methods of the invention, including the 
methods that have been summarized in the paragraphs immediately above. 
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In additional significant embodiments of the display means displaying a representation of the 
extent of relatedness, the cells in at least one class obtained as described in the NUCLEIC ACID ANALYSIS 
section above. 

The present invention additionally provides a display means displaying a representation of the 
5 correlation between a first class of cells and a second class of cells. The cells in each class are chosen 
from among cells of a given cell type, cells from a given tissue, and cells from a given organ, as 
described above. The correlation reflects differences between the first class and the second class in the 
presence of a pair of nucleotide subsequences, each pair consisting of a first subsequence and a second 
subsequence and the nucleotide length separating the first and second subsequences of the pair, and in a 
10 quantitation of the extent to which each pair having the determined length is present in the cells. 

In an advantageous embodiment of this display means, the correlation is related to a set of 
orthonormal eigenvectors. The elements of the basis set upon which the eigenvectors are constructed 
reflect particular biochemical or physiological pathways correlated between the cells of the two classes 
Each of these eigenvectors is associated with an eigenvalue that is an integer greater than zero. After 
1 5 defining an upper limit of the eigenvalues to be used, the coefficients of the basis set elements in each 
eigenvector whose eigenvalue is less than or equal to this upper limit reflect the contribution of the 
corresponding pathway to the biochemical or physiological differences correlated between the cells of 
the first class and the cells of the second class- 

In an additional advantageous embodiment of the display means displaying a representation of 
20 the correlation, the representation is a cluster diagram or a dendrogram, and includes a tree structure 
reflecting the relatedness of the pathways involved in the biochemical or physiological response to a 
difference between cells of the two classes. In obtaining this representation, a correlation matrix is 
calculated that provides a distance determination in which the distance reflects the amplitude of a 
difference vector. This vector is a difference between two vectors each of which reflects information 
25 obtained for the response of one of the two classes to the difference between the classes. The branches 
of the tree structure reflect the difference vectors and the branches are ramified from nodes. 

In important embodiments of the display means displaying a representation of the correlation, 
the representation is obtained employing the methods of the invention, including the methods that have 
been summarized in the paragraphs immediately above, 

30 In important embodiments of the representation of the correlation, the representation is obtained 

employing the methods of the invention, including the methods that have been summarized in the 
paragraphs immediately above. 



WO 00/15851 .23 - PCT/US99/21525 

In additional significant embodiments of the display means displaying a representation of the 
correlation, the cells in at least one class obtained as described in the NUCLEIC ACID ANALYSIS section 
above. 

Other Aspects 

In addition to providing representations of cells, the techniques described here are also useful for 
providing representations of nucleic acid fragments or genes. The starting point for the analysis is the 
matrix I sd described previously, where s labels the sample (or group of samples or distinct types of cells) 
and d labels a particular measurement of the expression level of a particular gene in that class. Rather 
than generating representations based on the rows of I, each representing a different sample or group of 
samples, it is possible to generate representations based on the columns of I, each representing a different 
nucleic acid. Hierarchical and geometrical representations of nucleic acids, based on their relative 
abundance across a series of cells, can be used to infer genes that are co-expressed and are likely to have 
related biological function. 

Other Embodiments 

The data matrix of intensities I can be described more generally as a representation in which 
each row corresponds to a particular biological sample or group of samples, and each column 
corresponds to a particular nucleic acid molecule or class of molecules whose quantities are measured in 
each of the biological states . 

In addition to the differential-display methods described to provide measurements of nucleic acid 
quantities, other methods for obtaining measurements of the nucleic acids present in a cell are available. 
These include restriction fragment length polymorphism, amplification fragment length polymorphism, 
EST sequencing, serial analysis of gene expression, hybridization to oligonucleotide probes, and other 
methods known in the art. Other methods, such as quantification by TaqMan or Northern blots, are also 
used. All of these methods generate data sets that can be analyzed according to the methods described 
here. The measurements I sd for each biological state and nucleic acid can correspond to absolute 
concentrations, concentrations relative to a standard (either ratio or numeric difference), or other 
convenient measures. 

The methods of the invention includes analysis of populations ranging from 5, 10, 25, 50, 100, 
1000, 10,000 or 100,000 or more members. 
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EXAMPLE 

Male Sprague-Dawley rats (Harlan Sprague Dawley, Inc., Indianapolis, Indiana) of 10-14 weeks 
of age were gavage-fed and dosed once a day for three days with the following drugs, dissolved in sterile 
water, at the following levels: 

5 phenobarbitol 3.81 mg/kg/day 

gabapentin 34.29 mg/kg/day 

vigabatrin 150 mg/kg/day 

paraldehyde 77.08 mg/kg/day. 

These dosages correspond to the EDI 00 (the upper limit of the effective dose for humans) 
10 adjusted for the difference in metabolic rate between rats and humans. Three rats were used for each 

drug treatment, and an additional three rats to match each drug were treated with sterile water to serve as 
a control. 

Rats were sacrificed 24 hours after the final dose and their brains were harvested. Collection of 
mRNA, synthesis of cDNA, and differential display protocols were carried out according to methods 
15 described in U. S. Patent No. 5,871,697 and Shimkets et al 1999 (Nature Biotechnology 1 7:798-803). 

The following steps were followed to analyze the differential display pattern: 

1 . The intensities I sr (x) for each of the three animals treated with the same drug were combined 
into a single average IJx), where the subscript a labels the drug. The standard deviation s ar (x) was also 
computed for the measurements from the individual animals treated with the drug. 

20 2. The averages I ar (x) and standard deviations s ar (x) for each drug were compared with the 

average I cr (x) and standard deviation s cr (x) for the sterile water control treatment, A difference at length 
x was marked if 

ABS(ln[I ar (x)/I cr (x)])>ln(L5) (5) 

and if the significance was smaller than 0.1 5 for a two-tailed t-test with 

25 t = [IJx) - IJx)] i [{ sjxf + s cr (x) 2 }/ 2 ] 1/2 (6) 

and infinite degrees of freedom. The difference intensities marked according to this procedure 
may then be inspected by eye and visually significant differences may be retained. 
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3. For each of the differences d, defined by a restriction enzyme pair r and a position x, the 
intensity I ar (x) = I atl was determined for each of the drug treatments, whether or not that particular 
treatment has a difference compared to the control. 

In this example, the final data matrix I ad has 8 rows: 1 row for each of the 4 drugs, and 1 row for 
5 each of 4 replicates of the water control data„ The matrix has as many columns as the number of 
differences detected in the differential display pattern. 

The Pearson correlation coefficient C ab between the 8 classes of samples (4 drugs, 4 water 
controls) was determined using methods provided in the Detailed Description of the Invention. If a data 
element for a particular difference was missing for a particular treatment, that difference did not 
1 0 contribute to the correlation coefficient. The correlations are shown in the Table 1 below, with the 
standard deviation within a drug shown as the diagonal elements. 

Table 1. 



15 





vigabatrir* <pheno barbito I 


gabapentin 


paraldehyde 


waterl 


water_2 


water_3 ; 


water_4 


:vigabatrin 


2613601! 


0.9982 


0.9913 


0.9728 : 


06548 


0.6673 


0 3674 


0 3786 


pheno barbital 


0 9982 263 4469 


03973: 


0 9325 


05678 


04724 


0.6569 


0 6601 


gabapentin 


0.9913. 


0.9973 


556 8114 


0 9922 


0.6177 


0 6057 


0 3400 


0.2704 


paraldehyde 


0.9728 


09325 


0 9922 


423.1916 ; 


0.6485 


05307 


0 5952 


0.5702 


:water_1 


0 6548 


0 5678 


0.6177 


0 6485 


59 4573 


0 9718 


09896 


03735 


water_2 


0.6673 


0.4724 


0,6057 


0.5307 


0.9718 


62 3568 


0 9960 


03886 


water. 3 


0.3674 


0 6569 


0 3400 


0 5952 


6,9896 


09960 


107.2630 


0.9975 


: water_4 


0 3786 


0,6601 


0.2704 


0.5702 


0.9735 


0 9886 


0 9978 


123.0743 


Next the pairwise Pearson 


distance was calculated 


as described previously. The 


distance \ 


is shown in Table 2 below. 
























Table 2. 












vigabatrin phenobarbitoi 


gabapentin 


paraldehyde 


waterl 


water_2 


water_3 


water_4 


vigabatrin 


0 0000 


0 0597 


0 1319 


0.2333 


0 8309 


0 8157 


1.1248 


11148 


phenobarbitoi 


0.0597 


0 0000 


0.0741 


0 3675 


0 9297 


10272 


0 8284 


0 8245 


gabapentin 


0.131$ 


0 074 1 


0 0000 


0 1245 


0.8744 


0.8880 


1.1489 


12080 


paraldehyde 


0 2333 


0 3675 


0 1245 : 


0.0000 


0.8384 


03688 


0.8997 


0.9271 


water_1 


0 830$ 


0 9297 


0.8744 


0.8384 : 


0.0000 


0.2376 


0.1446 


0,2302 


; water_2 


0 8157 


10272 


0.8B80 


0 9688 


0.2376 


0.0000 


0.0895 


0 150-3 


: water_3 


1.1248 


0 8284 


1 1489 


0.8997 


0 1446 


0 0895 


0 0000 


0.0663 


water_4 


1.1148 


0.8245 


1 2080 


09271 


0.2302 


0 1509 


0 0663 


0 0000 



The distances were then used as input to a nearest-neighbor clustering algorithm. The resulting 
clusters, using sterile H 2 0 as an outgroup, was shown in Fig. 3. The horizontal distances in Fig. 1 A were 
20 proportional to the pairwise Pearson distance between clusters. 

The correlation matrix C ab also served as the starting point for principal factor analysis. First, 
principal components were calculated using the inner product matrix from multidimensional scaling 



B = HCH 



(2) 
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where C is the correlation matrix and II is the centering matrix. The k 1h principal component is 
then the k th eigenvector of B normalized to unit length and ordered by decreasing eigenvalue A, k , and the 
k th principal factor was obtained by scaling the eigenvector by X k m . Projections of the treatments and 
controls onto principal factors are shown in Table 3 below. 

5 Table 3. 



facto r: 


1 


2; 


3 


4 


5 


6 


7 




eigenvalue: 


1.841 


0 347 


0 093 


0.036 


0 00? 


0 000 


-0036 


-038& 


vigabatrin 


-0 513 


-0/63 


-0.099 


0.083 


-0.004 


0 000 


0 000 


0 000 


phenobarbitol 


-0.397: 


0 317 


-0.132 


-0 038 


-0 0« 


0.000 


0 000 


0,000 


gabapentin 


-0 580 


-0.157 


0 028 


-0.094 


-0 004 


0,000 ; 


0 000 


0.000 


paraldehyde 


-0 404 


0.127 


o m 


0.057 


0034 


0 000 


0,000 


0.000 


water_1 


0 368 


-0 «9 


0 109 


0.0 T7 


•0.060 


0 000 


oooo 


0000 


water_2 


0 422 


-0.300 


-0.091 


-0.017 


0 03 & 


0000 


0.000 


0,000 


water_3 


0.542 


0.156: 


0 04& 


-0,090 


0.012 


0.000 


0 000 


0.000 


water_4 


0 561 


0 190 


-0 05& ; 


0 082 


0 001 


0000 


0 000 


0000 



The components are ordered from l (most infonnative) to 8 (least informative). The negative 
eigenvalues arise from the method used to account for missing data. If missing data had been handled in 
an alternate manner, for example if a missing element had been set to the average value or if the analysis 
10 were restricted to differences for which no data was missing, the eigenvalues would all be non-negative. 

In Fig. 4, the treatments are displayed by projection onto principal factors. Factor 1 
discriminates between drugs, where it has a negative value, and controls, where it has a positive value. 
Factor 2 discriminates between the drug treatments.. 

EQUIVALENTS 

15 From the foregoing detailed description of the specific embodiments of the invention, it should 

be apparent that unique methods for representing the extent of reiatedness between cells, cell lines, 
tissues, organs, or expressed sequences based on a genomic analysis of gene expression have been 
described. Although particular embodiments have been disclosed herein in detail, this has been done by 
way of example for purposes of illustration only, and is not intended to be limiting with respect to the 

20 scope of the appended claims which follow. In particular, it is contemplated by the inventor that various 
substitutions, alterations, and modifications may be made to the invention without departing from the 
spirit and scope of the invention as defined by the claims. For instance, the choice of source material, 
subsequences used, or software algorithm used is believed to be a matter of routine for a person of 
ordinary skill in the art with knowledge of the embodiments described herein* 
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CLAIMS 

I claim; 

1 . A method for generating a representation of the extent of reiatedness between at least two 
classes of cells, wherein the cells in each class are chosen from the group consisting of cells of a given 
cell type, cells from a given tissue, and cells from a given organ, the method comprising the steps of 

a) defining a plurality of pairs of nucleotide subsequences, each pair consisting of a first 
subsequence and a second subsequence; 

b) in the nucleic acid of each class of cells determining the presence of a fragment with the first 
subsequence at one end and the second subsequence at another end and having a length separated by the 
first and second subsequences, and a quantitation of the extent to which each fragment is present; and 

c) determining the extent of reiatedness reflecting similarities or differences in the presence and 
quantitation of the fragments among the classes. 

2. The method described in claim 1 wherein the determining of the presence and quantitation of 
the fragments described in step b) is carried out by a process comprising the steps of: 

i) digesting samples of the nucleic acid from the cells of each class with a plurality of specific 
pairs of restriction endonucleases, each sample being treated by one pair, one nuclease of the pair 
targeting the first subsequence and the second nuclease of the pair targeting the second subsequence, 
each digestion providing specific restriction fragments, hybridizing double stranded adapter DNA 
molecules to the fragments, each adapter DNA molecule comprising (a) a shorter strand having no 5* 
terminal phosphate and consisting of a first and second portion, said first portion being at the 5 5 end and 
being complementary to the overhang produced by one of the restriction endonucleases of the pair, and 
(b) a longer strand having a V end complementary to the second portion of the shorter strand, and 
ligating the longer strands to the fragments to produce ligated fragments, wherein each Iigated fragment 
is capable of generating an output signal; 

ii) generating output signals from each ligated fragment for each of the pairs of restriction 
endonucleases, each output signal characterizing (a) the subsequences of the pairs of restriction 
endonucleases (b) the length between the two subsequences corresponding to the two restriction 
endonucleases employed in each pair of nucleases, and (c) the quantitation of the fragment 
corresponding to the pair and the length; and 
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iii) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (a) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (b) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 

thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
cells in each class contains fragments having the specific subsequence pairs and the nucleotide length 
between the pairs. 

3. The method described in claim 1 wherein the determining of the presence of the fragments and 
the quantitation of the fragments, described in step b) is carried out by a process comprising the steps of: 

i) for each pair of nucleotide subsequences providing a pair of oligonucleotide primers, 
consisting of a first primer and a second primer, wherein the first primer is complementary to the first 
subsequence and the second primer is complementary to the second subsequence; 

ii) amplifying the nucleotide sequence between the first subsequence and the second 
subsequence using the oligonucleotide primers to prime the amplification, providing an amplicon 
characterized by the subsequence pair, a length between the two subsequences corresponding to the two 
primers employed in each pair and a quantitation of the extent to which each amplicon is present; and 

iii) generating output signals for each amplicon, each output signal characterizing (a) the 
subsequences of the pairs of primers, (b) the length, and (c) the quantitation; and 

iv) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (a) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (a) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 
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thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
cells in each class contains the specific subsequence pairs and the nucleotide length between the pairs. 

4. The method described in claim 1 wherein the extent of relatedness in step c) is provided by 
calculating a distance wherein the distance reflects the amplitude of a difference vector that is a 
difference between a first vector which reflects information derived from the quantitation for each 
subsequence pair obtained for the first class and a second vector which reflects information derived from 
the quantitation for each subsequence pair obtained for the second class, wherein different elements of 
each vector relate to data obtained using different pairs. 

5. The method described in claim 1 wherein the extent of relatedness in step c) is provided by 
generating a tree structure reflecting the relatedness between any two classes, wherein the branches of 
the tree structure reflect the difference vectors and the branches are ramified from nodes. 

6. The method described in claim 1 wherein the cells in at least one class are cancer cells. 

7. The method described in claim 1 wherein the cells in at least one class have been contacted 
with a putative pharmaceutical agent. 

8. A method for generating a representation of the correlation between a plurality of classes of 
cells wherein the cells in each class are chosen from the group consisting of cells of a given cell type, 
cells from a given tissue, and cells from a given organ, the correlation reflecting a change in the nature 
and amount of nucleic acids present in the classes, the method comprising the steps of: 

a) defining a plurality of pairs of nucleotide subsequences, each pair consisting of a first 
subsequence and a second subsequence; 

b) in the nucleic acid of each class of cells determining the presence of a fragment with the first 
subsequence at one end and the second subsequence at another end and having a length separated by the 
first and second subsequences, and a quantitation of the extent to which each fragment is present, thereby 
defining a difference between the classes; 

c) evaluating the correlation between the cells of the classes; and 

d) preparing a representation of the correlation. 
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9. The method described in claim 8 wherein the determining of the presence and quantitation of 
the fragments described in step b) is carried out by a process comprising the steps of: 

i) digesting samples of the nucleic acid from the cells of each class with a plurality of specific 
pairs of restriction endonucleases, each sample being treated by one pair, one nuclease of the pair 
targeting the first subsequence and the second nuclease of the pair targeting the second subsequence, 
each digestion providing specific restriction fragments, hybridizing double stranded adapter UNA 
molecules to the fragments, each adapter DNA molecule comprising (a) a shorter strand having no 5' 
terminal phosphate and consisting of a first and second portion, said first portion being at the 5' end and 
being complementary to the overhang produced by one of the restriction endonucleases of the pair, and 
(b) a longer strand having a 3' end complementary to the second portion of the shorter strand, and 
ligating the longer strands to the fragments to produce I i gated fragments, wherein each iigated fragment 
is capable of generating an output signal; 

ii) generating output signals from each Iigated fragment for each of the pairs of restriction 
endonucleases, each output signal characterizing (a) the subsequences of the pairs of restriction 
endonucleases (b) the length between the two subsequences corresponding to the two restriction 
endonucleases employed in each pair of nucleases, and (c) the quantitation of the fragment 
corresponding to the pair and the length; and 

iii) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (a) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (b) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 

thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
ceils in each class contains fragments having the specific subsequence pairs and the nucleotide length 
between the pairs, 

1 0. The method described in claim 8 wherein the determining of the presence of the fragments 
and the quantitation of the fragments, described in step b) is carried out by a process comprising the steps 
of: 
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i) for each pair of nucleotide subsequences providing a pair of oligonucleotide primers, 
consisting of a first primer and a second primer, wherein the first primer is complementary to the first 
subsequence and the second primer is complementary to the second subsequence; 

ii) amplifying the nucleotide sequence between the first subsequence and the second 
subsequence using the oligonucleotide primers to prime the amplification, providing an amplicon 
characterized by the subsequence pair, a length between the two subsequences corresponding to the two 
primers employed in each pair and a quantitation of the extent to which each amplicon is present; and 

iii) generating output signals for each amplicon, each output signal characterizing (a) the 
subsequences of the pairs of primers, (b) the length, and (c) the quantitation; and 

iv) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (a) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (b) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 

thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
cells in each class contains the specific subsequence pairs and the nucleotide length between the pairs. 

1 1 . The method described in claim 8 wherein the correlation in step d) is related to a set of 
orthonormal eigenvectors, the elements of the basis set upon which the eigenvectors are constructed 
reflecting particular biochemical or physiological pathways correlated between the cells of the two 
classes, each eigenvector having an eigenvalue that is an integer greater than zero, the coefficients of the 
basis set elements in each eigenvector whose eigenvalue is less than or equal to a particular integer that 
is an upper limit of the eigenvalues used reflecting the contribution of the corresponding pathway to the 
biochemical or physiological differences correlated between the cells of the first class and the cells of the 
second class. 
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12. The method described in claim 8 wherein the representation is a cluster diagram or a 
dendrogram, includes a tree structure reflecting the relatedness of the pathways involved in the 
biochemical or physiological response to a difference between cells of the two classes, wherein a 
correlation matrix provides a distance determination wherein the distance reflects the amplitude of a 
difference vector that is a difference between two vectors each of which reflects information obtained for 
the response of one of the two classes to the difference, and wherein the branches of the tree structure 
reflect the difference vectors and the branches are ramified from nodes. 

13. The method described in claim 8 wherein the cells in at least one class are cancer cells. 

14. The method described in claim 8 wherein the cells in at least one class have been contacted 
with a putative pharmaceutical agent, and the method comprises the steps of: 

a) treating the cells of at least one class with an amount of the agent sufficient to effect a change 
in the state of those cells or with an amount of the agent less than or equal to a predetermined upper limit 
of dosing concentration; 

b) defining a plurality of pairs of nucleotide subsequences, each pair consisting of a first 
subsequence and a second subsequence; 

c) in the nucleic acid of each class of cells determining the presence of a fragment with the first 
subsequence at one end and the second subsequence at another end and having a length separated by the 
first and second subsequences, and a quantitation of the extent to which each fragment is present, thereby 
defining an effect of the agent; 

d) evaluating the correlation between the effect of the agent on the cells of the first class and the 
effect of the agent on the cells of another class; and 

e) preparing a representation of the correlation. 

15. A display means displaying a representation of the extent of relatedness between at least two 
classes of cells, wherein the cells in each class are chosen from the group consisting of cells of a given 
cell type, cells from a given tissue, and cells from a given organ, the extent of relatedness reflecting, in 
the nucleic acids of the classes of cells, similarities or differences in the presence of pairs of nucleotide 
subsequences, each pair consisting of a first subsequence and a second subsequence, a nucleotide length 
separating the first and second subsequences of the pair and a quantitation of the extent to which each 
pair having the determined length is in the classes of cells„ 
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1 6. The display means described in claim 15 wherein the extent of relatedness is related to a 
distance wherein the distance reflects the amplitude of a difference vector that is a difference between a 
first vector which reflects information derived from the quantitation for each subsequence pair obtained 
for the first class and a second vector which reflects information derived from the quantitation for each 
subsequence pair obtained for the second class, wherein different elements of each vector relate to data 
obtained using different pairs. 

1 7. The display means described in claim 1 5 wherein the representation includes a tree structure 
reflecting the relatedness between any two classes, and wherein the branches of the tree structure reflect 
the difference vectors and the branches are ramified from nodes. 

1 8. The display means described in claim 15 wherein the extent of relatedness is obtained by a 
process comprising the steps of 

a) defining a plurality of pairs of nucleotide subsequences, each pair consisting of a first 
subsequence and a second subsequence; 

b) in the nucleic acid of each class of cells determining the presence of a fragment with the first 
subsequence at one end and the second subsequence at another end and having a length separated by the 
first and second subsequences, and a quantitation of the extent to which each fragment is present; and 

c) determining the extent of relatedness reflecting similarities or differences in the presence and 
quantitation of the fragments among the classes 

19. The display means described in claim 18 wherein the determining of the presence and 
quantitation of the fragments described in step b) is carried out by a process comprising the steps of: 

i) digesting samples of the nucleic acid from the cells of each class with a plurality of specific 
pairs of restriction endonuc leases, each sample being treated by one pair, one nuclease of the pair 
targeting the first subsequence and the second nuclease of the pair targeting the second subsequence, 
each digestion providing specific restriction fragments, hybridizing double stranded adapter DNA 
molecules to the fragments, each adapter DNA molecule comprising (a) a shorter strand having no 5' 
terminal phosphate and consisting of a first and second portion, said first portion being at the 5' end and 
being complementary to the overhang produced by one of the restriction endonucleases of the pair, and 
(b) a longer strand having a 3 ? end complementary to the second portion of the shorter strand, and 
ligating the longer strands to the fragments to produce iigated fragments, wherein each ligated fragment 
is capable of generating an output signal; 
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ii) generating output signals from each Iigated fragment for each of the pairs of restriction 
endonucleases, each output signal characterizing (a) the subsequences of the pairs of restriction 
endonucleases (b) the length between the two subsequences corresponding to the two restriction 
endonucleases employed in each pair of nucleases, and (c) the quantitation of the fragment 
corresponding to the pair and the length; and 

iii) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (a) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (b) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 

thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
cells in each class contains fragments having the specific subsequence pairs and the nucleotide length 
between the pairs.. 

20. The display means described in claim 18 wherein the determining of the presence of the 
fragments and the quantitation of the fragments, described in step b) is carried out by a process 
comprising the steps of: 

i) for each pair of nucleotide subsequences providing a pair of oligonucleotide primers, 
consisting of a first primer and a second primer, wherein the first primer is complementary to the first 
subsequence and the second primer is complementary to the second subsequence; 

ii) amplifying the nucleotide sequence between the first subsequence and the second 
subsequence using the oligonucleotide primers to prime the amplification, providing an amplicon 
characterized by the subsequence pair, a length between the two subsequences corresponding to 

the two primers employed in each pair and a quantitation of the extent to which each amplicon is 
present; and 

iii) generating output signals for each amplicon, each output signal characterizing (a) the 
subsequences of the pairs of primers, (b) the length, and (c) the quantitation; and 
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iv) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (a) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (a) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 

thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
cells in each class contains the specific subsequence pairs and the nucleotide length between the pairs.. 

21 . The display means described in claim 15 wherein the cells in at least one class are cancer 

cells. 

22„ The display means described in claim 15 wherein the cells in at least one class have been 
contacted with a putative pharmaceutical agent. 

23. A display means displaying a representation of the correlation between a plurality of classes 
of cells, wherein the cells in each class are chosen from the group consisting of cells of a given cell type, 
cells from a given tissue, and cells from a given organ, the correlation reflecting, in the nucleic acids of 
the classes of cells, differences in the presence of a pair of nucleotide subsequences, each pair consisting 
of a first subsequence and a second subsequence and the nucleotide length separating the first and second 
subsequences of the pair, and a quantitation of the extent to which each pair having the determined 
length is present in the cells, between the classes. 

24. The display means described in claim 23 wherein the correlation is related to a set of 
orthonormal eigenvectors, the elements of the basis set upon which the eigenvectors are constructed 
reflecting particular biochemical or physiological pathways correlated between the cells of the two 
classes, each eigenvector having an eigenvalue that is an integer greater than zero, the coefficients of the 
basis set elements in each eigenvector whose eigenvalue is less than a particular integer that is chosen to 
be an upper limit of the eigenvalues reflecting the contribution of the corresponding pathway to the 
biochemical or physiological differences correlated between the cells of the first class and the cells of the 
second class. 
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25. The display means described in claim 23 wherein the representation is a cluster diagram or a 
dendrogram and includes a tree structure reflecting the relatedness of the pathways involved in the 
biochemical or physiological difference between cells of the two classes, wherein a correlation matrix 
provides a distance determination wherein the distance reflects the amplitude of a difference vector that 
is a difference between two vectors each of which reflects information obtained for the difference 
between the classes, and wherein the branches of the tree structure reflect the difference vectors and the 
branches are ramified from nodes. 

26. The display means described in claim 23 wherein the correlation is obtained by a method 
comprising the steps of 

a) defining a plurality of pairs of nucleotide subsequences, each pair consisting of a first 
subsequence and a second subsequence; 

b) ) in the nucleic acid of each class of cells determining the presence of a fragment with the first 
subsequence at one end and the second subsequence at another end and having a length separated by the 
first and second subsequences, and a quantitation of the extent to which each fragment is present, thereby 
defining a difference between classes; 

c evaluating the correlation between the cells of one class and the cells of a second class based 
on the difference between them; and 

d) ) preparing a representation of the correlation . 

27. The display means described in claim 23 wherein the determining of the presence and 
quantitation of the fragments described in step b) is carried out by a process comprising the steps of: 

i) digesting samples of the nucleic acid from the cells of each class with a plurality of specific 
pairs of restriction endonucleases, each sample being treated by one pair, one nuclease of the pair 
targeting the first subsequence and the second nuclease of the pair targeting the second subsequence, 
each digestion providing specific restriction fragments, hybridizing double stranded adapter DNA 
molecules to the fragments, each adapter DNA molecule comprising (a) a shorter strand having no 5' 
terminal phosphate and consisting of a first and second portion, said first portion being at the 5' end and 
being complementary to the overhang produced by one of the restriction endonucleases of the pair, and 
(b) a longer strand having a 3' end complementary to the second portion of the shorter strand, and 
li gating the longer strands to the fragments to produce ligated fragments, wherein each li gated fragment 
is capable of generating an output signal; 
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ii) generating output signals from each ligated fragment for each of the pairs of restriction 
endonucleases, each output signal characterizing (a) the subsequences of the pairs of restriction 
endonucleases (b) the length between the two subsequences corresponding to the two restriction 
endonucleases employed in each pair of nucleases, and (c) the quantitation of the fragment 
corresponding to the pair and the length; and 

iii) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (a) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (b) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 

thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
cells in each class contains fragments having the specific subsequence pairs and the nucleotide length 
between the pairs. 

28, The display means described in claim 23 wherein the determining of the presence of the 
fragments and the quantitation of the fragments, described in step b) is carried out by a process 
comprising the steps of: 

i) for each pair of nucleotide subsequences providing a pair of oligonucleotide primers, 
consisting of a first primer and a second primer, wherein the first primer is complementary to the first 
subsequence and the second primer is complementary to the second subsequence; 

ii) amplifying the nucleotide sequence between the first subsequence and the second 
subsequence using the oligonucleotide primers to prime the amplification, providing an amplicon 
characterized by the subsequence pair, a length between the two subsequences corresponding to the two 
primers employed in each pair and a quantitation of the extent to which each amplicon is present; and 

iii) generating output signals for each amplicon, each output signal characterizing (a) the 
subsequences of the pairs of primers, (b) the length, and (c) the quantitation; and 
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iv) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (a) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (a) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 

thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
cells in each class contains the specific subsequence pairs and the nucleotide length between the pairs. 

29. The display means described in claim 23 wherein the cells in at least one class are cancer 

cells. 

30. The display means described in claim 23 wherein the cells in at least one class have been 
contacted with a putative pharmaceutical agent, and the correlation is obtained by method comprising the 
steps of 

a) contacting the cells of at least one class with an amount of the agent sufficient to effect a 
change in the state of those cells or with an amount of the agent less than or equal to a predetermined 
upper limit of dosing concentration; 

b) defining a plurality of pairs of nucleotide subsequences, each pair consisting of a first 
subsequence and a second subsequence; 

c) in the nucleic acid of each class of cells determining the presence of a fragment with the first 
subsequence at one end and the second subsequence at another end and having a length separated by the 
first and second subsequences, and a quantitation of the extent to which each fragment is present, thereby 
defining an effect of the agent; 

d) evaluating the correlation between the effect of the agent between the cells of at least one 
class contacted with the agent and the cells of another class; and 

e) preparing a representation of the correlation. 
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3 1 . A representation of the extent of relatedness between at least two classes of cells, wherein the 
cells in each class are chosen from the group consisting of cells of a given cell type, cells from a given 
tissue, and cells from a given organ, the extent of relatedness reflecting, in the nucleic acids of the 
classes of cells, similarities or differences in the presence of pairs of nucleotide subsequences, each pair 
consisting of a first subsequence and a second subsequence, a nucleotide length separating the first and 
second subsequences of the pair and a quantitation of the extent to which each pair having the 
determined length is in the classes of cells* 

32.. The representation described in claim 3 1 wherein the extent of relatedness is related to a 
distance wherein the distance reflects the amplitude of a difference vector that is a difference between a 
first vector which reflects information derived from the quantitation for each subsequence pair obtained 
for the first class and a second vector which reflects information derived from the quantitation for each 
subsequence pair obtained for the second class, wherein different elements of each vector relate to data 
obtained using different pairs. 

33. The representation described in claim 3 1 wherein the representation includes a tree structure 
reflecting the relatedness between any two classes, and wherein the branches of the tree structure reflect 
the difference vectors and the branches are ramified from nodes. 

34. The representation described in claim 31 wherein the extent of relatedness is obtained by a 
process comprising the steps of 

a) defining a plurality of pairs of nucleotide subsequences, each pair consisting of a first 
subsequence and a second subsequence; 

b) in the nucleic acid of each class of cells determining the presence of a fragment with the first 
subsequence at one end and the second subsequence at another end and having a length separated by the 
first and second subsequences, and a quantitation of the extent to which each fragment is present; and 

c) determining the extent of relatedness reflecting similarities or differences in the presence and 
quantitation of the fragments among the classes 
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35. The representation described in claim 34 wherein the determining of the presence and 
quantitation of the fragments described in step b) is carried out by a process comprising the steps of: 

i) digesting samples of the nucleic acid from the cells of each class with a plurality of specific 
pairs of restriction endonucleases, each sample being treated by one pair, one nuclease of the pair 
targeting the first subsequence and the second nuclease of the pair targeting the second subsequence, 
each digestion providing specific restriction fragments, hybridizing double stranded adapter DNA 
molecules to the fragments, each adapter DNA molecule comprising (a) a shorter strand having no 5' 
terminal phosphate and consisting of a first and second portion, said first portion being at the 5' end and 
being complementary to the overhang produced by one of the restriction endonucleases of the pair, and 
(b) a longer strand having a V end complementary to the second portion of the shorter strand, and 

li gating the longer strands to the fragments to produce li gated fragments, wherein each li gated fragment 
is capable of generating an output signal; 

ii) generating output signals from each ligated fragment for each of the pairs of restriction 
endonucleases, each output signal characterizing (a) the subsequences of the pairs of restriction 
endonucleases (b) the length between the two subsequences corresponding to the two restriction 
endonucleases employed in each pair of nucleases, and (c) the quantitation of the fragment 
corresponding to the pair and the length; and 

iii) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (a) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (b) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 

thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
cells in each class contains fragments having the specific subsequence pairs and the nucleotide length 
between the pairs. 

36. The representation described in claim 34 wherein the determining of the presence of the 
fragments and the quantitation of the fragments, described in step b) is carried out by a process 
comprising the steps of: 
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i) for each pair of nucleotide subsequences providing a pair of oligonucleotide primers, 
consisting of a first primer and a second primer, wherein the first primer is complementary to the first 
subsequence and the second primer is complementary to the second subsequence; 

ii) amplifying the nucleotide sequence between the first subsequence and the second 
subsequence using the oligonucleotide primers to prime the amplification, providing an amplicon 
characterized by the subsequence pair, a length between the two subsequences corresponding to 

the two primers employed in each pair and a quantitation of the extent to which each amplicon is 
present; and 

iii) generating output signals for each amplicon, each output signal characterizing (a) the 
subsequences of the pairs of primers, (b) the length, and (c) the quantitation; and 

iv) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (a) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (a) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 

thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
cells in each class contains the specific subsequence pairs and the nucleotide length between the pairs. 

37. The representation described in claim 3 1 wherein the cells in at least one class are cancer 

cells. 

38. Tlfie representation described in claim 31 wherein the cells in a class have been contacted 
with a putative pharmaceutical agent. 
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38. A representation of the correlation between a plurality of classes of cells, wherein the cells in 
each cla^s are chosen from the group consisting of cells of a given cell type, cells from a given tissue, 
and cells from a given organ, the correlation reflecting, in the nucleic acids of the classes of cells, 
differences in the presence of a pair of nucleotide subsequences, each pair consisting of a first 
subsequence and a second subsequence and the nucleotide length separating the first and second 
subsequences of the pair, and a quantitation of the extent to which each pair having the determined 
length is present in the cells, between the classes. 

39. The representation described in claim 38 wherein the correlation is related to a set of 
orthonormal eigenvectors, the elements of the basis set upon which the eigenvectors are constructed 
reflecting particular biochemical or physiological pathways correlated between the cells of the two 
classes, each eigenvector having an eigenvalue that is an integer greater than zero, the coefficients of the 
basis set elements in each eigenvector whose eigenvalue is less than a particular integer that is chosen to 
be an upper limit of the eigenvalues reflecting the contribution of the corresponding pathway to the 
biochemical or physiological differences correlated between the cells of the first class and the cells of the 
second class. 

40. The representation described in claim 38 wherein the representation is a cluster diagram or a 
dendrogram and includes a tree structure reflecting the reJatedness of the pathways involved in the 
biochemical or physiological differences between cells of the two classes, wherein a correlation matrix 
provides a distance determination wherein the distance reflects the amplitude of a difference vector that 
is a difference between two vectors each of which reflects information obtained from one of the classes, 
and wherein the branches of the tree structure reflect the difference vectors and the branches are 
ramified from nodes . 

4L The representation described in claim 38 wherein the correlation is obtained by a method 
comprising the steps of 

a) defining a plurality of pairs of nucleotide subsequences, each pair consisting of a first 
subsequence and a second subsequence; 

b) in the nucleic acid of each class of cells determining the presence of a fragment with the first 
subsequence at one end and the second subsequence at another end and having a length separated by the 
first and second subsequences, and a quantitation of the extent to which each fragment is present, thereby 
defining a difference between classes; 
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c) evaluating the correlation between the cells of one class and the cells of a second class based 
on the difference between them; and 

d) preparing a representation of the correlation. 

42. The representation described in claim 41 wherein the determining of the presence and 
quantitation of the fragments described in step b) is carried out by a process comprising the steps of: 

i) digesting samples of the nucleic acid from the cells of each class with a plurality of specific 
pairs of restriction endonucleases, each sample being treated by one pair, one nuclease of the pair 
targeting the first subsequence and the second nuclease of the pair targeting the second subsequence, 
each digestion providing specific restriction fragments, hybridizing double stranded adapter DNA 
molecules to the fragments, each adapter DNA molecule comprising (a) a shorter strand having no 5' 
terminal phosphate and consisting of a first and second portion, said first port ion being at the 5' end and 
being complementary to the overhang produced by one of the restriction endonucleases of the pair, and 
(b) a longer strand having a 3 ? end complementary to the second portion of the shorter strand, and 
ligating the longer strands to the fragments to produce li gated fragments, wherein each li gated fragment 
is capable of generating an output signal; 

ii) generating output signals from each ligated fragment for each of the pairs of restriction 
endonucleases, each output signal characterizing (a) the subsequences of the pairs of restriction 
endonucleases (b) the length between the two subsequences corresponding to the two restriction 
endonucleases employed in each pair of nucleases, and (c) the quantitation of the fragment 
corresponding to the pair and the length; and 

Hi) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (a) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (b) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 

thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
cells in each class contains fragments having the specific subsequence pairs and the nucleotide length 
between the pairs. 
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43. The representation described in claim 41 wherein the determining of the presence of the 
fragments and the quantitation of the fragments, described in step c) is carried out by a process 
comprising the steps of: 

i) for each pair of nucleotide subsequences providing a pair of oligonucleotide primers, 
consisting of a first primer and a second primer, wherein the first primer is complementary to the first 
subsequence and the second primer is complementary to the second subsequence; 

ii) amplifying the nucleotide sequence between the first subsequence and the second 
subsequence using the oligonucleotide primers to prime the amplification, providing an amplicon 
characterized by the subsequence pair, a length between the two subsequences corresponding to the two 
primers employed in each pair and a quantitation of the extent to which each amplicon is present; and 

iii) generating output signals for each amplicon, each output signal characterizing (a) the 
subsequences of the pairs of primers, (b) the length, and (c) the quantitation; and 

iv) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (a) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (a) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 

thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
cells in each class contains the specific subsequence pairs and the nucleotide length between the pairs, 

44. The representation described in claim 38 wherein the cells in at least one class are cancer 

cells. 

45. The representation described in claim 38 wherein the cells in at least one class have been 
contacted with a putative pharmaceutical agent, and the correlation is obtained by a method comprising 
the steps of 

a) contacting the cells of at least one class with an amount of the agent sufficient to effect a 
change in the state of those cells or with an amount of the agent less than or equal to a predetermined 
upper limit of dosing concentration; 
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b) defining a plurality of pairs of nucleotide subsequences, each pair consisting of a first 
subsequence and a second subsequence; 

c) in the nucleic acid of each class of cells determining the presence of a fragment with the first 
subsequence at one end and the second subsequence at another end and having a length separated by the 
first and second subsequences, and a quantitation of the extent to which each fragment is present, thereby 
defining an effect of the agent; 

d) evaluating the correlation between the effect of the agent between the cells of at least one 
class contacted with the agent and the cells of another class; and 

e) preparing a representation of the correlation. 

46. A method for generating a geometrical representation between a plurality of classes of cells 
wherein the cells in each class are chosen from the group consisting of cells of a given cell type, cells 
from a given tissue, and cells from a given organ, the representation reflecting a change in the nature and 
amount of nucleic acids present in the classes, the method comprising the steps of: 

a) in the nucleic acid of each class of cells, assessing the presence and amount of a nucleic acid 
fragment thereby defining a difference between the classes; 

b) carrying out a geometrical analysis based on the differences between the cells of the classes; 

and 

c) preparing a representation of the results of the analysis. 

47. The method described in claim 46 wherein the geometrical representation is a result obtained 
by a principal component analysis or a principal factor analysis. 

48. The method described in claim 46 wherein assessing the presence and amount of a nucleic 
acid fragment described in step a) is carried out by a process comprising the steps of: 

i) probing the nucleic acid of each class with a set of oligonucleotide probes specific for the 
fragment; and 

ii) determining the extent to which each probe binds the nucleic acid; 

thereby providing an assessment of the presence and amount of the nucleic acid fragment in the 

class. 

49. The method described in claim 46 wherein assessing the presence and amount of a nucleic 
acid fragment described in step a) is carried out by a process comprising the steps of:: 
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i) defining a plurality of pairs of nucleotide subsequences, each pair consisting of a first 
subsequence and a second subsequence, and 

ii) in the nucleic acid of each class of cells determining the presence of a fragment with the first 
subsequence at one end and the second subsequence at another end and having a length separated by the 
first and second subsequences, and a quantitation of the extent to which each fragment is present, thereby 
defining the difference between the classes. 

50. The method described in claim 49 wherein assessing the presence and quantity of a nucleic 
acid fragment described in step ii) is carried out by a process comprising the steps of: 

(a) digesting samples of the nucleic acid from the cells of each class with a plurality of specific 
pairs of restriction endonucleases, each sample being treated by one pair, one nuclease of the pair 
targeting the first subsequence and the second nuclease of the pair targeting the second subsequence, 
each digestion providing specific restriction fragments, hybridizing double stranded adapter DNA 
molecules to the fragments, each adapter DNA molecule comprising (1) a shorter strand having no 5' 
terminal phosphate and consisting of a first and second portion, said first portion being at the 5' end and 
being complementary to the overhang produced by one of the restriction endonucleases of the pair, and 
(2) a longer strand having a 3 ' end complementary to the second portion of the shorter strand, and 
ligating the longer strands to the fragments to produce ligated fragments, wherein each ligated fragment 
is capable of generating an output signal; 

(b) generating output signals from each ligated fragment for each of the pairs of restriction 
endonucleases, each output signal characterizing ( 1 ) the subsequences of the pairs of restriction 
endonucleases (2) the length between the two subsequences corresponding to the two restriction 
endonucleases employed in each pair of nucleases, and (3) the quantitation of the fragment 
corresponding to the pair and the length; and 

(c) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (1) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (2) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 



WO 00/15851 - 47 PCT/US99/21525 

thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
cells in each class contains fragments having the specific subsequence pairs and the nucleotide length 
between the pairs. 

51. The method described in claim 49 wherein assessing the presence and quantity of a nucleic 
acid fragment described in step ii) is carried out by a process comprising the steps of: 

(a) for each pair of nucleotide subsequences providing a pair of oligonucleotide primers, 
consisting of a first primer and a second primer, wherein the first primer is complementary to the first 
subsequence and the second primer is complementary to the second subsequence; 

(b) amplifying the nucleotide sequence between the first subsequence and the second 
subsequence using the oligonucleotide primers to prime the amplification, providing an amplicon 
characterized by the subsequence pair, a length between the two subsequences corresponding to the two 
primers employed in each pair and a quantitation of the extent to which each amplicon is present; and 

(c) generating output signals for each amplicon, each output signal characterizing (1) the 
subsequences of the pairs of primers, (2) the length, and (3) the quantitation; and 

(d) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (1) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (2) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 

thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
cells in each class contains the specific subsequence pairs and the nucleotide length between the pairs. 

52. The method described in claim 46 wherein the results of the geometrical analysis are chosen 
from the group consisting of eigenvalues, eigenvectors, and principal factors . 
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53. The method described in claim 46 wherein the results of analysis in step c) are related to a set 
of orthonormal eigenvectors, the elements of the basis set upon which the eigenvectors are constructed 
reflecting particular biochemical, physiological or pharmacological components correlated between the 
cells of the two classes, each eigenvector having an eigenvalue, the coefficients of the basis set elements 
in each eigenvector reflecting the contribution of the corresponding biochemical, physiological or 
pharmacological components to the differences between the cells of the first class and the cells of the 
second class. 

54. The method described in claim 46 wherein the cells in at least one class are cancer cells. 

55. The method described in claim 46 wherein the cells in at least one class are contacted with a 
putative pharmaceutical agent, and the method comprises the steps of: 

a) treating the cells of at least one class with an amount of the agent sufficient to effect a change 
in the state of those cells or with an amount of the agent less than or equal to a predetermined upper limit 
of dosing concentration; 

b) defining a plurality of pairs of nucleotide subsequences, each pair consisting of a first 
subsequence and a second subsequence; 

c) in the nucleic acid of each class of cells determining the presence of a fragment with the first 
subsequence at one end and the second subsequence at another end and having a length separated by the 
first and second subsequences, and a quantitation of the extent to which each fragment is present, thereby 
defining an effect of the agent; 

d) conducting a principal component analysis between the effect of the agent on the cells of the 
first class and the cells of another class; and 

e) preparing a representation of the results of the analysis . 

56. A display means displaying a geometrical representation between a plurality of classes of 
cells wherein the cells in each class are chosen from the group consisting of cells of a given cell type, 
cells from a given tissue, and cells from a given organ, the principal component analysis reflecting a 
change in the nature and amount of nucleic acids present in the classes, wherein the representation is 
obtained by a method comprising the steps of: 

a) in the nucleic acid of each class of cells, assessing the presence and amount of a nucleic acid 
fragment thereby defining a difference between the classes; 

b) carrying out a principal component analysis based on the differences between the ceils of the 
first class and the cells of the second class; and 
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c) preparing the representation of the results of the analysis, 

57. The display means described in claim 56 wherein the geometrical representation is a result 
obtained by a principal component analysis or a principal factor analysis. 

58. The display means described in claim 56 wherein assessing the presence and amount of a 

! 

nucleic acid fragment described in step a) comprises the steps of: 

i) probing the nucleic acid of each class with a set of oligonucleotide probes specific for the 
fragment; and 

ii) determining the extent to which each probe binds the nucleic acid; 

thereby providing an assessment of the presence and amount of the nucleic acid fragment in the 

class. 

59. The display means described in claim 56 wherein assessing the presence and amount of a 
nucleic acid fragment described in step a) is carried out by a process comprising the steps of: 

i) defining a plurality of pairs of nucleotide subsequences, each pair consisting of a first 
subsequence and a second subsequence; and 

ii) in the nucleic acid of each class of cells determining the presence of a fragment with the first 
subsequence at one end and the second subsequence at another end and having a length separated by the 
first and second subsequences, and a quantitation of the extent to which each fragment is present, thereby 
defining the difference between the classes. 

60. The display means described in claim 59 wherein determining the presence and quantity of a 
nucleic acid fragment described in step ii) is carried out by a process comprising the steps of: 

(a) digesting samples of the nucleic acid from the cells of each class with a plurality of specific 
pairs of restriction endonucleases, each sample being treated by one pair, one nuclease of the pair 
targeting the first subsequence and the second nuclease of the pair targeting the second subsequence, 
each digestion providing specific restriction fragments, hybridizing double stranded adapter DNA 
molecules to the fragments, each adapter DNA molecule comprising (1) a shorter strand having no 5' 
terminal phosphate and consisting of a first and second portion, said first portion being at the 5' end and 
being complementary to the overhang produced by one of the restriction endonucleases of the pair, and 
(2) a longer strand having a 3' end complementary to the second portion of the shorter strand, and 
ligating the longer strands to the fragments to produce ligated fragments, wherein each ligated fragment 
is capable of generating an output signal; 
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(b) generating output signals from each ligated fragment for each of the pairs of restriction 
endonuc leases, each output signal characterizing ( 1 ) the subsequences of the pairs of restriction 
endonucleases (2) the length between the two subsequences corresponding to the two restriction 
endonucleases employed in each pair of nucleases, and (3) the quantitation of the fragment 
corresponding to the pair and the length; and 

(c) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (1) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (2) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 

thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
cells in each class contains fragments having the specific subsequence pairs and the nucleotide length 
between the pairs. 

61 . The display means described in claim 59 wherein assessing the presence and quantity of a 
nucleic acid fragment described in step ii) is carried out by a process comprising the steps of: 

(a) for each pair of nucleotide subsequences providing a pair of oligonucleotide primers, 
consisting of a first primer and a second primer, wherein the first primer is complementary to the first 
subsequence and the second primer is complementary to the second subsequence; 

(b) amplifying the nucleotide sequence between the first subsequence and the second 
subsequence using the oligonucleotide primers to prime the amplification, providing an amplicon 
characterized by the subsequence pair, a length between the two subsequences corresponding to the two 
primers employed in each pair and a quantitation of the extent to which each amplicon is present; and 

(c) generating output signals for each amplicon, each output signal characterizing (1) the 
subsequences of the pairs of primers, (2) the length, and (3) the quantitation; and 
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(d) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (1) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (2) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 

thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
cells in each class contains the specific subsequence pairs and the nucleotide length between the pairs. 

62. The display means described in claim 56 wherein the results of the analysis are chosen from 
the group consisting of eigenvalues, eigenvectors, and principal factors. 

63. The display means described in claim 56 wherein the results of the analysis in step c) are 
related to a set of orthonormai eigenvectors, the elements of the basis set upon which the eigenvectors 
are constructed reflecting particular biochemical, physiological or pharmacological components 
correlated between the cells of the two classes, each eigenvector having an eigenvalue, the coefficients of 
the basis set elements in each eigenvector reflecting the contribution of the corresponding biochemical, 
physiological or pharmacological components to the differences between the cells of the first class and 
the cells of the second class. 

64. The display means described in claim 56 wherein the cells in at least one class are cancer 

cells. 

65. The display means described in claim 56 wherein the cells in at least one class have been 
contacted with a putative pharmaceutical agent, and the representation is obtained by a method 
comprising the steps of: 

a) treating the cells of at least one class with an amount of the agent sufficient to effect a change 
in the state of those cells or with an amount of the agent less than or equal to a predetermined upper limit 
of dosing concentration; 

b) defining a plurality of pairs of nucleotide subsequences, each pair consisting of a first 
subsequence and a second subsequence; 
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c) in the nucleic acid of each class of cells determining the presence of a fragment with the first 
subsequence at one end and the second subsequence at another end and having a length separated by the 
first and second subsequences, and a quantitation of the extent to which each fragment is present, thereby 
defining an effect of the agent; 

d) conducting a principal component analysis between the effect of the agent on the cells of the 
first class and the cells of another class; and 

e) prepar ing the representation of the results of the analysis. 

66. A geometrical representation between a plurality of classes of cells wherein the cells in each 
class are chosen from the group consisting of cells of a given cell type, cells from a given tissue, and 
cells from a given organ, the principal component analysis reflecting a change in the nature and amount 
of nucleic acids present in the classes, the representation obtained by a method comprising the steps of; 

a) in the nucleic acid of each class of cells, assessing the presence and amount of a nucleic acid 
fragment thereby defining a difference between the classes; 

b) carrying out a principal component analysis based on the differences between the cells of the 
first class and the cells of the second class; and 

c) preparing the representation of the results of the analysis. 

67. The representation described in claim 66 wherein the geometrical representation is a result 
obtained by a principal component analysis or a principal factor analysis. 

68 . The representation described in claim 66 wherein assessing the presence and amount of a 
nucleic acid fragment described in step a) comprises the steps of: 

i) probing the nucleic acid of each class with a set of oligonucleotide probes specific for the 
fragment; and 

ii) determining the extent to which each probe binds the nucleic acid; 

thereby providing an assessment of the presence and amount of the nucleic acid fragment in the 

class. 
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69. The representation described in claim 66 wherein assessing the presence and amount of a 
nucleic acid fragment described in step a) is carried out by a process comprising the steps of: 

i) defining a plurality of pairs of nucleotide subsequences, each pair consisting of a first 
subsequence and a second subsequence; and 

ii) in the nucleic acid of each class of cells determining the presence of a fragment with the first 
subsequence at one end and the second subsequence at another end and having a length separated by the 
first and second subsequences, and a quantitation of the extent to which each fragment is present, thereby 
defining the difference between the classes. 

70. The representation described in claim 69 wherein determining the presence and quantity of a 
nucleic acid fragment described in step ii) is carried out by a process comprising the steps of: 

(a) digesting samples of the nucleic acid from the cells of each class with a plurality of specific 
pairs of restriction endonucleases, each sample being treated by one pair, one nuclease of the pair 
targeting the first subsequence and the second nuclease of the pair targeting the second subsequence, 
each digestion providing specific restriction fragments, hybridizing double stranded adapter DNA 
molecules to the fragments, each adapter DNA molecule comprising (1) a shorter strand having no 5' 
terminal phosphate and consisting of a first and second portion, said first portion being at the 5' end and 
being complementary to the overhang produced by one of the restriction endonucleases of the pair, and 
(2) a longer strand having a 3' end complementary to the second portion of the shorter strand, and 
ligating the longer strands to the fragments to produce ligated fragments, wherein each ligated fragment 
is capable of generating an output signal; 

(b) generating output signals from each ligated fragment for each of the pairs of restriction 
endonucleases, each output signal characterizing ( 1 ) the subsequences of the pairs of restriction 
endonucleases (2) the length between the two subsequences corresponding to the two restriction 
endonucleases employed in each pair of nucleases, and (3) the quantitation of the fragment 
corresponding to the pair and the length; and 
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(c) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (1) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (2) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 

thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
cells in each class contains fragments having the specific subsequence pairs and the nucleotide length 
between the pairs. 

7L The representation described in claim 69 wherein assessing the presence and quantity of a 
nucleic acid fragment described in step ii) is carried out by a process comprising the steps of: 

(a) for each pair of nucleotide subsequences providing a pair of oligonucleotide primers, 
consisting of a first primer and a second primer, wherein the first primer is complementary to the first 
subsequence and the second primer is complementary to the second subsequence; 

(b) amplifying the nucleotide sequence between the first subsequence and the second 
subsequence using the oligonucleotide primers to prime the amplification, providing an amplicon 
characterized by the subsequence pair, a length between the two subsequences corresponding to the two 
primers employed in each pair and a quantitation of the extent to which each amplicon is present; and 

(c) generating output signals for each amplicon, each output signal characterizing (1 ) the 
subsequences of the pairs of primers, (2) the length, and (3) the quantitation; and 

(d) optionally searching a nucleotide sequence database to determine sequences that are 
predicted to produce or the absence of any sequences that are predicted to produce the one or more 
output signals produced by the nucleic acid from the cells of each class, the database comprising a 
plurality of known nucleotide sequences of nucleic acids that may be present in the cells of each class, a 
sequence from the database being predicted to produce the one or more output signals when the sequence 
from the database has both (1) the same length between occurrences of target nucleotide subsequences as 
is represented by the one or more output signals, and (2) the same target nucleotide subsequence as are 
represented by said one or more output signals, or target nucleotide subsequences that are members of 
the same stets of target nucleotide subsequences represented by the one or more output signals, 
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thereby providing a quantitative measure of the extent to which the nucleic acid present in the 
cells in each class contains the specific subsequence pairs and the nucleotide length between the pairs. 

12, The representation described in claim 66 wherein the results of the analysis are chosen from 
the group consisting of eigenvalues, eigenvectors, and principal factors. 

73. The representation described in claim 66 wherein the results of the analysis in step c) are 
related to a set of orthonormal eigenvectors, the elements of the basis set upon which the eigenvectors 
are constructed reflecting particular biochemical, physiological or pharmacological components 
correlated between the cells of the two classes, each eigenvector having an eigenvalue, the coefficients of 
the basis set elements in each eigenvector reflecting the contribution of the corresponding biochemical, 
physiological or pharmacological components to the differences between the ceils of the first class and 
the cells of the second class. 

74. The representation described in claim 66 wherein the celts in at least one class are cancer 

cells„ 

75. The representation described in claim 66 wherein the cells in at least one class have been 
contacted with a putative pharmaceutical agent, and the representation is obtained by a method 
comprising the steps of: 

a) treating the cells of at least one class with an amount of the agent sufficient to effect a change 
in the state of those cells or with an amount of the agent less than or equal to a predetermined upper limit 
of dosing concentration; 

b) defining a plurality of pairs of nucleotide subsequences, each pair consisting of a first 
subsequence and a second subsequence; 

c) in the nucleic acid of each class of cells determining the presence of a fragment with the first 
subsequence at one end and the second subsequence at another end and having a length separated by the 
first and second subsequences, and a quantitation of the extent to which each fragment is present, thereby 
defining an effect of the agent; 

d) conducting a principal component analysis between the effect of the agent on the cells of the 
first class and the cells of another class; and 

e) preparing the representation of the results of the analysis. 
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76. A method for classifying a plurality of classes of cells or components thereof hierarchically 
comprising the steps of 

a) measuring relative differences in the quantity of a nucleic acid present in each class of cells to 
provide measurements of differential nucleic acid display; 

b) converting the measurements into distances between the classes of cells in a vector space; and 

c) preparing a hierarchical classification amongst the classes based on the vector distances. 

77. The method of claim 76 wherein the classification is performed on classes of cells, wherein 
the cells in a class may be cells of a given cell type, cells from a given tissue, and cells from a given 
organ, cells exhibiting a particular pathological state, or cells which have been contacted with a putative 
pharmaceutical agent. 

78. The method of claim 76 wherein the classification is performed on a component of the cells 
in the classes, wherein the component comprises a gene, a nucleic acid, or a fragment thereof 

79. The method of claim 76 wherein the measuring is carr ied out by a procedure chosen from 
the group consisting of differential display of nucleic acid fragments, probing for the presence of a 
nucleic acid using an oligonucleotide probe, sequences obtained from expressed sequence tags (ESTs), 
assessing restriction fragment length polymorphisms, and assessing amplification fragment length 
polymorphisms 

80. The method of claim 76 wherein the preparation of the hierarchical classification is carried 
out by a procedure chosen from the group consisting of principal component analysis of a correlation 
matrix, principal factor analysis of a correlation matrix, principal component analysis of a centered inner 
product matrix, and principal factor analysis of a centered inner product matrix, 

8 1 . The method of claim 80 further comprising the step of obtaining a distance metric between 
the classes from a reduced dimensionality geometrical representation. 

82. A display means displaying the results of the classification obtained by a method described 
in any one of claims 76-8 L 

83 „ A method for representing a plurality of classes of ceils or components thereof geometrically 
comprising the steps of 

a) measuring relative differences in the quantity of a nucleic acid present in each class of cells to 
provide measurements of differential nucleic acid display; and 
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b) preparing a geometrical representation amongst the classes based on the measurement of the 
differential display. 

84. The method of claim 83 wherein the classification is performed on classes of cells, wherein 
the cells in a class may be cells of a given cell type, cells from a given tissue, and cells from a given 
organ, cells exhibiting a particular pathological state, or cells which have been contacted with a putative 
pharmaceutical agent. 

85. The method of claim 83 wherein the classification is performed on a component of the cells 
in the classes, wherein the component comprises a gene, a nucleic acid, or a fragment thereof. 

86. The method of claim 83 wherein the measuring is carried out by a procedure chosen from 
the group consisting of differential display of nucleic acid fragments, probing for the presence of a 
nucleic acid using an oligonucleotide probe, sequences obtained from expressed sequence tags (ESTs), 
assessing restriction fragment length polymorphisms, and assessing amplification fragment length 
polymorphisms 

87„ The method of claim 83 wherein the preparation of the hierarchical classification is carried 
out by a procedure chosen from the group consisting of principal component analysis of a correlation 
matrix, principal factor analysis of a correlation matrix, principal component analysis of a centered inner 
product matrix, and principal factor analysis of a centered inner product matrix. 

88. The method of claim 87 further comprising the step of obtaining a distance metric between 
the classes from a reduced dimensionality geometrical representation. 

89. A display means displaying the results of the geometrical representation obtained by a 
method described in any one of claims 83-88. 

90. A method of presenting the hierarchical relatedness of two or more members of a 
population, the method comprising: 

providing a data set of each member in the population; 

generating a hierarchical classification of said data set; and 

displaying said classification, thereby presenting the hierarchical relatedness of the members of 
the population. 

91. The method of claim 90, wherein said population is a population of cells. 

92.. The method of claim 90, wherein said population is a population of nucleic acid sequences. 
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93. The method of claim 90, wherein said population is a population of polypeptide sequences. 

94. The method of claim 90, wherein said hierarchical classification of any two or more 
members of the population is calculated using a distance method in combination with an algorithm* 

95 . The method of claim 94, wherein said distance method is a Pearson correlation distance, 
Euclidean distance, Manhattan distance, Mahalanobis distance, a pairwise Pearson distance, or a 
Spearman distance, 

96. The method of claim 95, wherein said algorithm is single linkage, average linkage, or 
complete linkage 

97. The method of claim 90, wherein said data set is the product of an analysis of said members 
of the population that is selected from the group consisting of differential display, serial analysis of gene 
expression , expression tagged sequence analysis, restriction fragment length polymorphism, amplified 
fragment length polymorphism, or Northern blot hybridization analysis. 

98. A method of presenting the geometrical relatedness of two or more members of a 
population, the method comprising: 

providing a data set of each member in the population; 

generating a geometrical classification of said data set; and 

displaying said classification, thereby presenting the geometrical relatedness of the members of 
the population. 

99. The method of claim 98, wherein said population is a population of cells. 

100. The method of claim 98, wherein said population is a population of nucleic acid sequences. 

101. The method of claim 98, wherein said population is a population of polypeptide sequences. 

102. The method of claim 98, wherein said geometrical classification is generated by analyzing 
a matrix using an algorithm. 

103. The method of claim 102, wherein said matrix includes a correlation matrix. 

104. The method of claim 103, wherein said correlation matrix includes a Pearson correlation 
matrix, a Spearman correlation matrix, or a pairwise Pearson correlation matrix. 

105. The method of claim 102, wherein said matrix includes a centered inner product distance 

matrix. 
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106. The method of claim 105, wherein the inner product distance matrix is determined using a 
distance calculated by hierarchical classification analysis „ 

107. The method of claim 102, wherein said algorithm includes principal component analysis . 

108. The method of claim 102, wherein said algorithm includes principal factor analysis. 

109. The method of claim 107, wherein said algorithm includes principal factor analysis. 

1 10„ The method of claim 102, wherein said geometrical classification is further analyzed using 
hierarchical classification. 

111. The method of claim 90, wherein said population includes 5, 10,25,50, 100, 1000, 10,000, 
100,000 or more members. 

112. The method of claim 98, wherein said population includes 5, 10, 25,50, 100, 1000, 10,000, 
1 00,000 or more members. 
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